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I. INTRODUCTION 



When analyzing or designing a control system, it is important to study the elTects 
on overall system performance of varying one or more parameters (mass, inertia, gain, 
resistance, etc). It is equally important to determine whether a desired dynamic behavior 
can be achieved with any set of values for the parameters ... if not, redesign is indicated. 

The analysis and synthesis of linear feedback control systems, or the compensation 
of same, can be realized by parameter plane methods. The parameter plane method, 
which works well for two variable parameters and which may be extended to three or 
more parameters, is purely algebraic, and the resulting plots are valuable aids to analysis. 

Furthermore, the performance index cost curves in terms of integral-square error 
and other cost functions can be superimposed on the parameter plane graphically. That 
is very useful for analysis and for designing a control system when its characteristic 
equation has two parameters. Since the coefficients of the characteristic polynomial are 
determined by the system parameters, it follows that some relationship exists between 
the value of any parameter and the value of the characteristic roots. The integral-square 
error as a performance index is also determined by the system parameters. 

Let us examine the choice of the integral-square error as a performance index. In 
many applications the performance of feedback control systems can be considered either 
satisfactory if the errors are below a certain limit, or unsatisfactory, if the errors are 
above a certain limit. The hypothesis that the configuration is fixed refers to the form 
of the system; it is assumed that one or more parameters can be adjusted. Given the 
input and ideal output signals as transient functions of time and given the configuration 
of the control system, find the values that the free parameters should have in order to 
minimize the integral-square error. 

One approach to solving this problem starts with the Fourier transforms of the input 
and ideal output signals. The integral-square error is then easily formulated by means 
of Parseval's theorem. If the transform of the error is a rational function of the complex 
frequency, the integral-square error can be easily evaluated by means of a table of in- 
tegral, which is in Appendix A. When the integral-square error is evaluated, it appears 
as a function of the free parameters. 

The object of this thesis, then, is to develop a user friendly interactive computer 
program to plot families of performance index cost curves. A secondary objective is to 
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demonstrate the use of these curves in the design of the system. The design procedure 
starts with the performance index cost curves superimposed on parameter plane curves. 
The result of the analysis is a choice of values for the parameters such that the system 
performance is the best available compromise between transient performance and cost. 
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II. BACKGROUND THEORY AND EQUATIONS 



A. FUNDAMENTAL CONCEPT 

Classical control theory, historically speaking, has been applied to the design prob- 
lem with the assistance of graphical presentation and trial and error methods. Such 
methods have been quite successful in the development of good control systems, but do 
not answer the question “Is this the best system possible.” The engineer and mathe- 
matician have alwavs looked for the “best” control svstem. Earlv methods used “stand- 
ard forms.” The concept of a “performance index” was invoked, also at an early point, 
for, if the chosen index could be satisfied, one could say that the resulting system was 
indeed the “best” as measured by the index. 



Let us take the performance index that is the integral-square value of the error be- 
tween the ideal output and the actual output. A simple feedback control system with 
adjustable parameters is assumed as presented in Figure 1. 




figure 1. Block diagram of feedback control system. 
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The performance index in terms of integral-square error is defined as: 




E 1 dt 



( 2 . 1 ) 



where PI ; performance index 

E ; input - actual output 

and the criterion was to design the system such that this index was minimized. For lin- 
ear systems this problem can be solved by applying Parseval's theorem [Ref. 1: p. 43] to 
obtain explicit relationships for system parameters. One should note that a system de- 
signed to minimize this performance index has a step response with considerable over- 
shoot and oscillation, and this may not be acceptable for some applications, but where 
such dynamic characteristics are suitable the design is straightforward. Other perform- 
ance indices which have been found useful are 



Both of these are nonlinear and do not have analytic solutions. Designs using such 
criteria were accomplished using simulation. The solution of an optimal control problem 
is the result of minimizing the chosen cost function and is the design which best satisfies 
the criteria defined by that cost function. 

B. PARSEVAL'S THEOREM 

Suppose we wish to know the area beneath the product of two functions .v,(r) and 
.r,(r) of the infinite interval of the time from minus infinity to plus infinity. Let this area 
be denoted by the symbol I (which stands for integral), thus 




( 2 . 2 ) 




(2.3) 




— oo 



(2.4) 
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Suppose, in addition, that the time functions are Fourier transformable and have trans- 
forms .V[(s) and .v,(s), respectively. We wish to express the integral I directly in terms 
of these transforms in order to obviate inverse transforming .v ( (s) and .v,(s) when these 
are given rather than the time functions themselves. To do this we observe that 



* 2 (0 = 



1 

2 7/ 



r+joo 



e SI .r 2 (s) ds 



-yoo 



Substitution of this value of x 2 {t) into Equation (2.4) yields 



(2.5) 




dt 



(2.6) 



Let us now interchange the order of the integrations so that we integrate with respect 
to time first 



1 f 4-y'oo p +oo 

I=TI- *2(0^5 e s, Xl (')dt (2.7) 

2n J —joo J -oo 

By the direct Fourier transform we know that 

p+oo 

*i(0- e~ s, x x {t)dt (2.8) 

^ — OO 

Thus the integral with respect to time on the right side of Equation (2.7) may be evalu- 
ated as 



p+oo 

e%(t)* = .v,(-s) (2.9) 

^ — OO 

This permits us to write Equation (2.7) as 

l f +7°° 

l= 2ri\ *i(“0*2(0<fc (2-10) 

^ *2 ^ —joo 

This is the desired expression for the integral I. 
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An important particular case of the above result occur when the two time functions 
involved are identical so that the integral I becomes 

r+oo 

1= x\t)dt (2.11) 

—GO 

We call this integral the integral-square value of the time function x(t). 

By Equation (2.10) we know this can be expressed as 



I = 




'+joo 

x( — s ) .v(s) ds 

J -Jvo 



( 2 . 12 ) 



This result is known as Parseval's theorem. This theorem represents a very convenient 
way of expressing the integral-square value of a time function in terms of its transform. 

C. MINIMIZATION PROCEDURE 

On the basis of Parseval's theorem we are now in a position to formulate the pro- 
cedure for minimizing the integral-square error for transient input signals by adjustment 
of the free parameters. 

This procedure involves four steps. The first step is to express the Fourier transform 
of the error as a function of the complex frequency s. This function will involve the free 
parameters of the system as unknown coefficients. The second step is to express the 
integral-square error I, in terms of the error transform E(s) by Parseval's theorem, i.e., 



= .f +70 ° £( — " )£W ds (2 - 13) 

At this stage, provided E(s) is a rational function, the integral-square error will appear 
in the form 



1 * 



j r+yoo 




c( -s) c(.s) 
d( — 5 ) d(s) 



(2.14) 



where c(s) and d(s ) are polynomials in s. The third step is to evaluate the integral. 
Fortunately, definite integrals of this form have been evaluated in terms of the coeffi- 
cients appearing in the polynomials. A table of these integrals along with an explanation 
of how they were evaluated is given in Appendix A. Upon carrying out the integration 
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we have the integral-square error expressed as a function of the free parameters 
p x through p L Symbolical!}' the integral-square error at this stage is expressed as 

l e (f\,P2, Pk ) 

The fourth step of the minimization procedure is to adjust the values of the free pa- 
rameters in such a way as to minimize the integral-square error. Formally this can be 
done in the standard manner of equating the partial derivatives of the integral-square 
error to zero and solving the resultant set of equations for the values of the parameters. 
In symbols this means solving the set of k equations of the form 



Unfortunately in many practical problems this formal procedure for finding the param- 
eter values that minimize the integral-square error leads to sets of non-linear equations 
in the parameters which are so complex that straightforward methods of solution do not 
exist. 

If the equations for the parameters obtained by setting the partial derivatives of the 
integral-square error with respect to the parameters equal to zero are too complex to 
be solved by analytical methods, we may resort to trial-and-error or successive approxi- 
mation techniques. However, if such methods are necessary to solve the equations for 
the parameters, it is probably best to go back to the expression for the integral-square 
error itself and plot this as a function of one of the parameters, holding the other pa- 
rameters fixed. By using a sufficient number of fixed values for these other parameters 
and making a sufficient number of plots, values can be found for all the parameters that 
define the minimum value of the integral-square error closely enough for engineering 
purposes. By working directly with the integral-square error expression, we avoid taking 
the partial derivative of this expression with respect to each of the free parameters. 
Taking these derivatives may be a very difficult and time-consuming chore. 
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III. PERFORMANCE INDEX CURVE EVALUATION 



A. BACKGROUND 

1. Overview 

The performance index program is an interactive personal computer program 
whenever possible. Before using this program, we have to generate a datafile containing 
a . [1 and PI. This file will enable us to plot families of performance index cost curves. 
There currently exists, at Naval Postgraduate School, a DSL(Dynamic Simulation Lan- 
guage) [Ref. 2] program residing on the school's IBM 370 mainframe. By using the DSL 
program, we can easily obtain the above mentioned data file, and then, down load this 
file from the mainframe computer to floppy diskette. This data, now on floppy diskette, 
is to be used, along with the computer program developed in this research paper, to plot 
performance index curves. With this program, the user can examine particular values of 
the performance index curves. Options, such as printer type and size of output, can be 
selected by the user. 

2. Software / Hardware 

The program language used in coding this program is Microsoft FORTRAX77 
V3.20. The Plotworks PLOT88 graphics library [Ref. 3] is used to generate output plots. 
The plotting of performance index curves conducted for this thesis were performed on 
an IBM-AT or IBM compatible 80286 machine. Due to the universality of FORTRAN 
coding, this program could be implemented on any machine capable of being pro- 
grammed in FORTRAN. The source code of generating curves of constant perform- 
ance index is listed in Appendix B. 

Before running the program. ANSI. SYS should be incorporated in the personal 
computer's CONFIG.SYS file prior to running the performance index program. 
Through the ANSI. SYS device driver, system calls to clear the screen and position the 
cursor are enabled. If not, screen readability suffers although the program is still fully 
functional. More will be said on this in the ANSI module section of this chapter. 



3. Two Dimensional Contouring 

The control grid is constructed as follows Figure 2. 
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Figure 2. The contour grid two-dimensional real array z from Ref.3. 
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where • = dimensional size of real array z. 

(nx.ny) = maximum index of area used. 

(xlow.ylow) = lower left corner of the grid area in the world units, 
(xhigh.yhigh) = upper right corner of the grid area in the world units, 
(nxsize.nysize) = grid interval in world units, 
dx = (xhigh-xlow),'(nx-l) 
dy = (yhigh-ylow) (ny-1) 



Note that only the control points in the range (i = 1 ,nx ; j = 1, ,ny) are used 

in the array z. By convention, the first index is referred to as the row index and the sec- 
ond index is referred to as the column index. 

The following sequence produces a two-dimensional contour map: 

• A grid constructed over the surface Figure 2. 

At each control point of the grid, the z ordinate is determined. Using Laplace or 
Spline interpolation or a combination of the two. curves of constant z are con- 
structed by PLOTSS graphics library. To do that, we can use the subroutine named 
ZGR1D which is the PLOTSS library, this subroutine interpolates a set of arbi- 
trarily spaced data points xp(k) ,yp(k), k = 1, npts onto a rectangular grid z{i,j) of 
arbitrary size. 

• The grid control points may be smoothed. The ZSMTH subroutine which is the 
PLOTSS library smooths the gridded control points by applying Laplacian 
smoothing. 

• Potions of the contour may be removed to make room for annotation. 

• The contours are drawn. Each set of 4 grid points can be subdivided further into 
4 or 16 subgrids for more precise contouring. In this program, the set of 4 grid 
points is subdivided into 16 subgrids. Variable contour line types and widths can 
be drawn. To make drawing more precisely, we can use the subroutine named 
ZCSEG which is the PLOTSS library. This subroutine converts the grid of control 
points into a set of contours and draws the contours. In addition, the capability 
of dividing each grid square into 1, 4, 16 subsections using cubic polynomial in- 
terpolation. This step may be needed if contours are close together or if the line 
quality of contour need improvement. 

The size of the contour grid is defined in the main program. For most applica- 
tions. a size of 50 x 50 control points is sufficient. The maximum number of control 
points is limited by the available space in the data segment of the program. A maximum 
of the 51 control levels can be drawn. 
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B. THE PROGRAM 



The performance index program given in Appendix B itself consists of main program 
and six subroutines which make up the performance index package. Figure 3 is a flow 
chart illustrating the organization. 




Figure 3. Performance Index Program Flow Chart. 



When the program is initiated, the following questions are presented: 
1. How many constant PI curves # = = = > 

• The user has to enter an integer value which is must be less than 51. 
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2. Enter PI values, lwgt, ldig # = = = > 
• Table 1 give us proper values. 



Table 1. PARAMETER TABLE 



parameter 


type 


value 


Description 


PI 


real 




Performance Index 


lwgt 


integer 


1 to 27 


the line type for each PI curve the PI curve 
is drawn with a solid line. See PLOTS8 
Manual. (Ref. 3: p. 32j 


ldit 


integer 


-I to -27 


the PI curves is drawn with a dashed line. 


ldit 


integer 


> 0 


the annotation control information for each 
PI. number of digits after decimal point. 


ldit 


integer 


-1 


omit decimal point. 


ldit 


integer 


< -2 


omit contour annotation for this level. 



3. Enter number of data file points = = = > 

• The user has to enter total number of data points. 

1. Subroutine SIZE 

The subroutine SIZE Figure 4 is displayed. This subroutine permits adjustment 
of plot size to both the screen and the printer. This option was utilized to scale plot 
outputs for inclusion in this thesis. 

*** This is the routine to Adjust the plot size *** 



Output plots are currently sized to fill the console screen. 

Would you like to adjust the plot size? 

Enter "y" or "n” ==> y 

A default plot factor size of 0.9 is used to fill the console screen. 
This plot size will be halved by entering a value of 0.45. 

It will be doubled with a factor entry of 1.8. 

Enter 'factor size ==> 0.7 



Figure 4. The plot size subroutine. 
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2. Subroutine LDFILE 

The subroutine LDflLE (Figure 5) is displayed on the screen. 



*** This is the routine to LOAD data from a file *** 



File name should not exceed 8 characters 
File extension should not exceed 3 characters 

What is the file name (fn) and extension (ext)? 

Enter in form fn.ext (e.g. PARAPLAN . INP ) ===> pil.dat 



Figure 5. The data load subroutine. 

Subroutine LDFILE permits program input of a pre-existing data file. The user 
is queried as to file name (not to exceed eight characters) and file extension (three char- 
acters or less). The code checks for existence of the file within the working subdirectory. 
If it does not exist, a message is returned stating such and offering the option of entering 
another file name again. Entry of file extension is optional. In fact, a file can be labeled 
with as little as one character or number. 

We note that the data should be formatted in order to read the data file suc- 
cessfully by using the LDFILE. A data file could be created aligning in the format 
identified by the subroutine LDFILE source code listed in Appendix B. 

3. Subroutine MONPRT 

The largest of the subroutines is subroutine MONPRT. A number of smaller 
subroutines are associated with it. Subroutine MONPRT allows selection of the graphics 
output device, be it the monitor or a wide variety of printers. Graphics output defaults 
to the monitor when the program is first executed. Graphics output can always be 
dumped from the screen to a printer, but there may be times when direct output to a 
printer is desired. All NTS printers are included in the printer output options. In addi- 
tion, one option permits the user to directly enter the values for IOPORT and MODEL, 
as outlined in the PLOT88 manual. [Ref. 3: p. 17) 
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The first menu displayed when subroutine MONPRT is called is the 



printer output menu Figure 6. 





[ PRINTER/OUTPUT MENU ! 


1 

1 


PRINTER 


NO. ! 


PRINTER 


i 

» 


1 

1 


1 


i 

i 


Epson FX-80 , All 


i 

i 


1 

1 


2 


i 

i 


Epson FX-100, All 


» 

i 


1 

1 


3 


i 

i 


Epson MX-100, All 


i 

i 


1 

1 


4 


i 

i 


Epson RX-80, All 


i 

t 


1 

1 


5 


i 

i 


Epson MX-80 & IBM Printer 


t 

i 


1 

1 


6 


i 

i 


HP 7470A Graphics Plotter 


i 

i 


1 

1 


7 


i 

i 


HP 7475A Graphics Plotter 


i 

i 


1 

1 


8 


i 

i 


HP 758xB Series Plotters 


i 

i 


1 

1 


9 


i 

i 


HP 2686A Laser Jet 


i 

i 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


1 

1 


10 


i 

i 


Graphics Monitor (default) 


i 

i 


1 

1 


11 


i 

i 


HARDWARE Interface Menu 


i 

i 


1 

1 


12 


i 

i 


Input PLOT88 Values for IOPORT and MODEL 


i 

i 


1 

1 


99 


i 

i 


RETURN to LOAD Data 


i 

i 




Enter 


integer 


number for selection = = = > 





Figure 6. The printer/output menu. 



Selection of a printer in this menu automatically directs graphics output to the most 
commonly associated output port, either parallel or serial, for that particular device. 
However, usually there are multiple parallel or serial ports to which a printer can be at- 
tached. The default printer output setting may not direct plots to the appropriate port. 
To provide maximum output flexibility, the printer/output menu allows selection of a 
specific output port through access to another menu. 
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| HARDWARE INTERFACE MENU | 


| SELECT NO. 


1 

1 


PORT 1 


1 1 


1 LPTl 


Printer Port J 


1 2 


LPT2 


Printer Port 


! 3 


| LPT3 


Printer Port | 


: 4 


1 COM1 


Serial Port J 


: 5 


! COM2 


Serial Port | 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


1 99 


1 RETURN to LOAD Data \ 



Enter integer number for selection ===> 



Figure 7. Tlie Hardware Interface Menu. 



The output port selection option is called the HARDWARE INTERFACE 
MENU Figure 7. This menu is part of a subordinate MONPRT subroutine titled 
PORT. The option allow selection of one of three parallel ports (LPTl - LPT3) or cither 
of two serial ports (COM1 - COM2). 

If a serial port is manually selected in the HARDWARE INTERFACE MENU, 
an associated data transfer (baud) rate must also be assigned. Immediately after se- 
lection of serial port options 4 or 5, the BAUD (data transfer) RATE MENU is dis- 
played Figure 8. Transfer of graphics data over a single line is usually time consuming, 
therefore the highest printer capable transfer rate available (9600 baud) is normally the 
best selection. A slower transfer rate would only be indicated if improper graphics out- 
put is generated. 



| BAUD (data transfer) RATE MENU ! 


1 SELECT NO. 


1 BAUD RATE ! 


! 1 


! 300 1 


1 2 


1 1200 ! 


1 3 


j 4800 1 


! 4 


1 9600 | 



Enter integer number for selection ===> 



Figure 8. Baud (data transfer) Rate Menu. 
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Bit by bit transfer of data can also have a check sum (parity) associated with it 
to provide an interval check of correct data transfer. Following selection of a data 
transfer rate, the PARITY MENU Figure 9 is displayed. The user has the option of 
choosing odd, even or no parity. 



! PARITY MENU ! 


! SELECT NO. i 


PARITY ■ J 


• 1 1 

1 A 1 


NO Parity J 


1 9 ■ 

1 ^ 1 


EVEN Parity | 


1 O 1 

1 ° 1 


ODD Parity [ 



Enter integer number for selection ===> 



Figure 9. Parity Menu. 



The PLOTS8 Software Library Reference Manual [Ref. 3: p. 18] contains a 
number of tables which provide common output settings for many more printers than 
can be selected in the PRINTER/OUTPUT menu. 

C. ANSI MODULE 

The ANSI module contains two subroutines. Subroutine CLR clears the monitor 
screen. It is used primarily as a precursor to displaying a menu or positioning queries 
or informational text. Subroutine CUP positions the cursor immediately after a request 
for response. This subroutine serves not only to give a pleasing screen display but also 
provide the user with an immediate reference as to where he is the program and if input 
is required of the user. Both routines can be found in the PLOT88 Reference Manual 
l Ref. 3]. The two subroutines which comprise the ANSI module make use of ANSI es- 
cape sequences. These sequences are a series of characters that can be used to define 
functions to MS-DOS. In addition to clearing the screen and positioning, moving, and 
saving the cursor's position, ANSI escape sequences can control screen graphics and 
reassign key definitions. 
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In order to utilize the escape sequences, the microprocessor must first be able to 
interpret the commands. It does this through a systems file entitled ANSI. SYS. 
ANSI. SYS must be defined as a device in the CONFIG.SYS file to enable proper inter- 
pretation of ANSI escape sequences (e.g., device= \ for \ ansi.sys where \ for \ is the 
path to ansi.sys). If ANSI.SYS is not resident in the CONFIG. SYS file, the perform- 
ance index program will not clear the screen or position the cursor when so directed. 
Instead, the indicated escape code sequence will be printed on the monitor whenever a 
call is made to CLR or CUP. Although program readability is degraded, the perform- 
ance index program will still function correctly. 
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IV. PARAMETER PLANE DESIGN -GRAPHICAL METHOD 



A. OVERVIEW 

The algebraic solutions have the disadvantage that a fixed value of zeta and omega 
must first be chosen to compute the alpha and beta terms. In some instances it is pos- 
sible to modify the remainder polynomial so as to ensure that the specified roots are 
dominant. However, it is not always possible to guarantee that roots placed at a speci- 
fied location can be made dominant. Thus, an exhaustive trial-and-error procedure may 
be required to achieve the best values for the various parameters. To avoid this trial- 
and-error analytical approach one can employ the graphical solution. 

We have two programs available. One permits us to get constant zeta curves and 
omega curves.The second program gives us a families of constant performance index 
curves in the parameter plane. We can incorporate the results of two programs by 
superimposing the plotted performance index cost curves over the parameter plane 
curves on the same sheet. 

To do that, we perform a preliminary procedure. This procedure involves four steps. 
The first step is to plot the parameter plane curves for the specific problem. This curve 
permits us to choose ranges of two parameters that are to be used in the analysis and 
design. The second step is to obtain a data file containing a , /? and performance index 
values by using the D.S.L program. We then copy this data file from the mainframe to 
a floppy diskette. The third step is to plot the families of performance index curves in 
the same range of parameter values on the two parameter plane. We then have two 
kinds of parameter plane curves. The fourth step of the analysis for a given specific 
problem is to duplicate these curves on the separate transparent sheets and then put 
them together on the same sheet. Recapitulation of the graphical analysis path is shown 
in Figure 10. 
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Figure M). Graphical analysis path. 



B. GRAF1I1CAL SOLUTIONS 

In this section, three dilfercnt kinds of control problems will be demonstrated. One 
of these is a cascade compensator system. The other one is a velocity and acceleration 
feedback compensator system. The above two problems will be evaluated with three 
kinds of performance index constraints, i.c., PI = f E 1 dt , V[—j\E\tdt and 
PI — f \ E\ dt. l hc third problem is a steering control for optimized autopilots. The 
controller parameters are adjusted to minimize the cost function. 

f \^ 2 e + d 2 )dt 

i J 0 



where <5 = rudder angle 
i//,= yaw error 
). - a weighting factor 
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1. Example 4-1 (cascade compensator) 

Figure 1 1 gives the block diagram of a high gain third order system which is to 
be damped with a cascade filter. (Ref. 4: p. 319) 




Figure 1 1. The third order cascade compensated system. 



Let 







The characteristic equation for this system is then: 

S 4 + (1 1 + /?)S 3 + ( 10+ lI/?)S 2 + (10/?+100cc)S + 100/? = 0 (4.1) 

Separating the coefficient terms. 

• Constant coefficient: 5 J + 11S 3 + 10S 2 

• ALPHA coefficient: 100S 

• BETA coefficient: S 3 + 11S 2 + 10S+100 

Enter each coefficient array value and desired zeta and omega values into the computer 
program, figure 12 provides a look at the combined zeta and omega graph spanning a 
frequency range from 1.0 to 30, C — 1.0 , 0.5 , 0.0 and co„ = 5.0 
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Figure 12. Combine zeta and omega curves for example 4.1 
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By inspection of the two parameter range, the alpha values vary from 0.0 to 140. 
and beta values also vary from 0.0 to 100. Now, within these ranges, simulation in 
D.S.L is used to evaluate the cost at a large number of points on the parameter plane. 
The D.S.L simulation program and its resulting data file containing a , /? and three dif- 
ferent criterion cost values are given in Appendix C. 

Once we have obtained a data file, we download this file from the mainframe 
computer to a lloppy diskette. To run the performance index program, type in the exe- 
cutable file name, COSTFUN. The IBM microcomputers in the NFS Controls Lab are 
set up under the file management system 1D1R. To execute the program under the this 
system, enter the appropriate subdirectory by using the up/down arrow keys to position 
the cursor adjacent to the subdirectory name, then depress the ENTER key. Use this 
same procedure to select COSTFUN.EXE and push ENTER to run the program. The 
preparation step for the plotting desired curves will appear. 

Enter proper values after each prompt as follows: 



1. How many constant PI curves # == 

2. Enter PI values ,lwgt,ldig # 1= 

2. Enter PI values ,lwgt,ldig # 2= 

2. Enter PI values ,lwgt,ldig # 3= 

2. Enter PI values ,lwgt,ldig # 4= 

2. Enter PI values ,lwgt,ldig # 5= 

2. Enter PI values ,lwgt,ldig # 6= 

3. Enter number of data file points 



>6 

=>0.27,1,2 
= > 0 . 2 , 1,1 
= >0 . 15,1,2 
=> 0 . 12 , 1,2 
=> 0 . 11 , 1,2 
= > 0.1 , 1,1 
===>86 



Next, adjust the plot size. 



*** This is the routine to Adjust the plot size *** 



Output plots are currently sized to fill the console screen. 

Would you like to adjust the plot size? 

Enter V or "n" ==> y 



A default plot factor size of 0.9 is used to fill the console screen. 
This plot size will be halved by entering a value of 0.45. 

It will be doubled with a factor entry of 1.8. 

Enter factor size ==> 0.7 



The last step is to enter a data file name. 
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*** This is the routine to LOAD data from a file *** 



File name should not exceed 8 characters 
File extension should not exceed 3 characters 



What is the file name { f n ) and extension (ext)? 

Enter in form fn.ext (e.g. PARAPLAN . INP ) ===> pil.dat 

Alter displaying the graph on the monitor, the program prompts you for a 
hardcopy of the graph. 

WOULD YOU LIKE HARDCOPY OF THESE PLOTS? (Y or N) 

y 



PRINTER/OUTPUT MENU 



PRINTER NO. 



1 

2 

3 

4 

5 

6 

7 

8 
9 

++++++++++++++ 



PRINTER 



! Epson FX-80, All 

! Epson FX-100, All ' 

1 Epson MX-100, All 

! Epson RX-80, All 

| Epson MX-80 & IBM Printer 

| HP 7470A Graphics Plotter 

| HP 7475A Graphics Plotter 

! HP 758xB Series Plotters 

! HP 2686A Laser Jet 

+++++++++++++++++++++++++++++++++++++++++++++++ 



10 

11 

12 

99 



Graphics Monitor (default) 

HARDWARE Interface Menu 

Input PLOT88 Values for IOPORT and MODEL 
RETURN to LOAD Data 



Enter integer number for selection ===> 5 

Reenter the data file name for the hard copy. Figure 13 is a performance index 
curve for constraint f E? dt . In like manner, the other two nonlinear criterion of per- 
formance index in terms of f\E\ ( dt and f\E\ dt can also be plotted by running the 
same program. Figures 14 and 15 are shown for desired constraint f\E\ t dt and 
f\E\ dt , respectively. 
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Figure 13. The performance index curves for f E 1 dt. 
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Figure 14. The performance index curves for f\E\ t dt. 
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5. The performance index curves for f\E\ dt. 
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Let us consider the two kinds of parameter plane curves, Figure 12 and 
Figure 13. By superimposing one over the other, the curves can produce the very in- 
teresting graph of Figure 16. In like manner, Figure 12 curves superimposed on the 
other cost contour are shown on Figure 17 and Figure 18. These graphs contain three 
kinds of constant families curves: £ , co n and cost curves. At this point, from Figure 16, 
we could select arbitrary operating points. The selected ALPHA BETA pairs are: 



A point: a = 25, /? = SO 



B point: a = 19.5, ft = 60 
C point: oc = 9.0, /? = 23 
D point: a = 100, /? = 70 

Now, go back to the characteristic equation Eq. (4.1) and substitute in this equation 
each ALPHA, BETA pair and then, simulate the time response for each one. The tab- 
ulated simulation results are shown in Figure 19 and Table 2 on page 32. From the 
simulation results, we find that operating point D (a = 25. /? = SO) had the fastest setting 
time, the lowest overshoot and the lowest cost index. The cost associated with this point 
are: 



E dt = 0. 1 1 



I E\ tdt — 0.12 
J | E| dt =0.23 

And step response is given in Figure 19. In this case, a cascade filter was designed with 
the zero equal to 3.2 ( = 80/25) and the pole equal to 80, using oc= P/Z and ft — P. 
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Figure 16. Combine zeta, omega and cost curves for f E? dt. 
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Figure 17. Combine zeta, omega and cost curves for f I E\ t dt. 
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Figure 18. Combine zeta, omega and cost curves for f | E \ dt. 
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Figure 19. Step response on example 4.1 
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Table 2. SIMULATION RESULTS. 



op- 

crating 

point 


parameters 


PI values 


K, 


t, 


a 


P 


JDJi 


f\E\ tdt 


f\E\ dt 


A 


25 


80 


0.27 


0.43 


0.55 


31 


2.3 


B 


19.5 


60 


0.27 


0.43 


0.53 


31 


2.2 


C 


9 


23 


0.3 


0.53 


0.6 


34 


1.9 


D 


100 


70 


0. 1 1 


0.12 


0.23 


29 


1.3 



The exact values, however, will be chosen to achieve a compromise between 
transient performance and cost, based on given system specifications. It is of interest 
to compare the performance index between the constraint f E? dt and the constraints 
f I E | t dt and f\ E\ dt at the same operating points. 



2. Example 4-2 (feedback compensator) 

This example deals with a plant incorporating both velocity (tachometer) and 
acceleration feedback. The block diagram of such a system is depicted in Figure 20. [Ref. 
4: p. 331[ 




Figure 20. Velocity and acceleration feedback compensated system. 
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The characteristic equation for this system is: 



S 3 + (3 + 50A2)S 2 + (2 + 50A1)S + 50 = 0 (4.2) 



Let 



a = K 1 



/? = A'2 

Using the same procedure as the previous example, Figure 21 provides the combined 
zeta and omega graph spanning a frequency range from 1.0 to 30. 

C - 0.45 

= 2-24 

This graph shows the two parameters range, the alpha values vary from 0.0 to 1.6 and 
beta values vary from 0.0 to 1.0. 

Now. we evaluate the cost at a large number of points on this range by using 
D.S.L. The D.S.L simulation program and its resultant data file containing a , /? and 
three different cost values are given in Appendix C. Once we have the data file, we can 
produce the family of performance index curves. These are shown in Figure 22, 
Figure 23 and Figure 24, respectively. 
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Figure 21. Constant zeta and omega curves for example 4.2 
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Figure 22. The performance index curves for / 1 ? dt. 



35 




Q 



o 

CJ 



6 



o 

CO 

■ 

o 



o 

CO 

■ 

n 



O 

=t< 

c? 



r-j 



o 



rd 

X 

VjL 



CE 



Figure 23. The performance index curves for f\E\ t dt. 
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Figure 24. 
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The performance index curves for f\E\ dt. 
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Let us consider the parameter plane curves. Figure 25, Figure 26 and Figure 27. These 
were produced by superimposing Figure 21 over each cost contour curve. Figure 22, 
Figure 23 and Figure 24. 

Figure 25 consists of three families of curves: zeta, omega and cost curves. We 
could select arbitrary operating points given by: 

A point: a = 0.45, /? = 0. 18 
B point: a = 0.60, /? = 0.40 
C point: a = 0.6, /? = 0.067 
D point: a = 0.8, /? = 0.22 
E point: a = 0.4, j3 = 0. 125 

The operating point C (a = 0.6. /? = 0.067) permits us the lowest cost index. The costs 
associated with this points are: 



E 2 dt = 0.47 
I E\ tdt =0.25 

J | £| dt = 0.68 

However, settling time was too slow, t s = 0.6 sec. In this system, operating point C 
(a = 0.6, /? = 0.067) determines the parameters K1 and K2, 

Kl = 0.6 

K2 = 0.067 
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Figure 25. Combine zeta, omega and cost curves for f t ? dt. 



39 



o 




®g 



Figure 26. Combine zeta, omega and cost curves for f\E\ t dt. 
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Figure 27. Combine zeta, omega and cost curves for f\E | dt. 
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Table 3. SIMULATION RESULTS. 



op- 

erating 

point 


parameters 


PI values 




t, 


a 


P 


jpd, 


f\E\ tdt 


/I El dt 


A 


0.45 


0.18 


0.6 


1.9 


1.1 


21 


4.1 


n 


0.6 


0.4 


0.76 


2.5 


1.4 


20 


5.9 


C 


0.6 


0.067 


0.47 


0.25 


0.68 


21 


6.0 


D 


0.8 


0.22 


0.62 


0.74 


0.95 


1 


2.1 


E 


0.4 


0.125 


0.6 


2.0 


1.1 


23 


4.5 



Table 3 and Figure 28 show simulation results for operating points A through E. The 
designer can then choose the operating point (a, /?) depending on his desired specifica- 
tions. 



3. Example 4-3 (optimization of controller) 

As a final engineering example, consider the steering control for a surface ship. 
The block diagram of such a system shown in Figure 29. 




Figure 29. Optimization of controller. 
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Given that the cost function is defined: 



J 




+ U 2 ) dt 



where <5 = rudder angle 
i j) e = yaw error 
/. = a weighting factor 

In this study, we used values given by R.E Reid [Ref. 5] for the SL-7 container 
ship. Reid uses a second order Nomoto model for the SL-7, 



<7(s) = 



K 

5(57" + 1) 



and also uses a controller 



G c (s) = 



Ail (s 71 + 1) 
{sT2 + 1) 



His results are given in Table 4. 



Table 4. REID RESULTS. 



ship 

speed 

Knots 


plant 


weight- 
ing factor 


controller 


K 


T 


K1 


T1 


T2 


16 


0.1 0S4 


90.36 


16.796 


0.4556 


89.51 


10.06 



Let us assume T1 and T2 are adjustable free parameters. 

r r 0.10S4 0.4556(571 + 1) 

c_ s(90.36s+ 1) * (sT2 + 1) 



Let 



a 



72’ ^ 72 

The characteristic equation for this system is: 



90.36S 3 + (1 + 90.36 p)S 2 + (/? +0.049a)S + 0.049£ = 0 



(4.3) 
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The parameter plane curves examined within the frequency range from 0.01 to 0.1 
C = 0.5, 1.0 window size; 0 < a < 35., 0 ^ /? <, 1 Using the parameter plane program, we 
get the praph of Figure 30. 



o 

o 







Figure 30. Constant zeta curves for example 4.3 
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In order to evaluate a family of cost curves, we have to obtain a data file on the same 
range of parameters, as Figure 30 by simulation. The D.S.L simulation program can 
be implemented with ease by rearranging the block diagram in Figure 31. 




Figure 31. Rearrangement of block diagram. 
The controller block becomes 



(* + 4 ) 

(s + /?) 



The D.S.L simulation program and data file for this example are given in APPENDIX 
C. Applying the previously discussed procedures for producing constant cost curves 
provides the graph of Figure 32. 
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From Figure 32. we could select the operating points given by: 

A point: a = 9, /? = 0. 1 
B point: <x = 20, /? = 0.2 
C point: a = 20. p = 0.4 

The operating point A (a = 9.0, /? = 0.1) provides the minimum cost. 
Recall that 





The two parameters values become 



T 1 = 90.0. T2 = 10.0 



Comparison of Reid's result and the simulation result show 'hat they are very closely 
matched. In addition, we simulated yaw response and rudder activity at the operating 
points A , B and C in Figure 33. Figure 34 compares the yaw response to course 
change for the ship at each operating point. Figure 35 compares the rudder activity. 
It is noted that operating point C is a faster controller but has a higher cost value. This 
cost is a higher rate of fuel consumption used to overcome the higher hydrodynamic 
drag forces. 
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Figure 33. Combine zeta and cost curves for -y- + <5 J ) dt. 
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Figure 3-1. The yaw response vs time. 
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Figure 35. The rudder activity vs time. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The graphical solution of the parameter plane technique provides the designer with 
a visual means of deducing how the dominant roots of the characteristic equation move 
about in the s-plane as two user-defined parameters defining compensator attributes are 
varied. This was illustrated in the example problems in Chapter 4. The graphical 
method is able to analyze easily, not only the linear constraint f E? di but also the non- 
linear constraint j\E\ t dt , f\E\ di and any kind of cost function using the perform- 
ance index program. 

Constant cost value contours can be plotted in much the same presentation as the 
parameter plane plot. This procedure is a useful design tool which permits simultaneous 
analysis of both dynamic response and cost. 

As has been demonstrated, rapid solution of a wide range of linear automatic con- 
trol systems problems can be achieved through use of the performance index program 
developed as part of this thesis. 



B. RECOMMENDATIONS 

Further work should be done to combine individual programs so that both param- 
eter plane and cost curves are produced by one consolidated program. 



APPENDIX A. TABULATED VALUES OF THE INTEGRAL FORM 



Table A. I gives the value of l n for values of rt from 1 to 9 where 

1 f + -/°° c(s)c( - 5 ) 

n 2 nj }_ joo aS d(s)d( s) 

and 

c(j) = c^,/ -1 + . . . + c 0 
d(s) = d n s + • • • + 



(A. 1-1) 



('Ll -2) 



(A. I — 3) 



Table A.l 



1 2 ^ ' 

_ Ci d n + c () d 2 

2 “ 2 dm ; 

J _ c 2^0^1 + ( C 1 ~ -CnC;Mp^3 + Cj^djdj 

3 “ 2d 0 d 3 ( -d 0 d 3 + <hd 2 ) ’ 

j _ C3 ( dfj d 3 + d () d ] d 2 ) + (c 2 — 2c ] c 3 )d () d ] d d + (C| — IcQC^dod^dj + c 0 ( —d^d^ + d 2 d 3 d d ) . 

2d 0 d^( —d§d 3 — d^d^ + d\d 2 d 3 ) 

I 5 = [_clm 0 + (c 3 - 2 c 2 c 4 )m, + (c 2 2 - 2c, c 3 + 2c 0 c 4 )m 2 + (c, 2 - 2c 0 c 2 )m 3 + c 0 2 m 4 ]; 
where 



,n ° = 



m, = —d 0 d 3 + d x d 2 . 
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f>h dqd$ 4 " d)d 4 , 



m 3 = -j- (d 2 m 2 - d A m x ), 
°o 



m A = — (d 2 m 3 — d A m 2 ), and 
d o 



A 5 = d Q {d } m A - d 3 m 3 + d 5 m 2 ). 



h = \_ c s m o + ( c l ~ 2c 3 c 5 )m, + (c 2 - 2c 2 c 4 + 2c,c s )w 2 + (c 2 - 2c, c 3 + 2c 0 c 4 )m 3 ] 

+ ~Ty [( c f ~ 2c 0 c 2 )m 4 c 0 2 m 5 ]; 



where 



m 0 = -j- {d A m x - d 2 m 2 + d 0 m 3 ), 
a 6 



o 2 

m, = —d 0 d ] d s + d 0 d 3 4- d ] d A — d l d 2 d 3 , 



2 

wij = d^d^d^ -h d l d 6 — 



9 3 

= -I - dfi d^d^ , 



m A = — {d 2 m 3 - d A tn 2 + d 6 m{), 
a o 



m 5 = (c/ 2 m 4 — ^/ 4 m 3 4- ^ 6 a? 2 2 ). and 



A 6 = ~ ^ 3 m 4 + d $ m 3 ). 



h = -^7 [ c 6 w o + i c l ~ 2c 4 c 6 )m, + (c 4 - 2c 3 c 5 + 2c 2 c 6 )w 2 + (c 3 - 2c 2 c 4 + 2c,c 5 - 2c 0 c 6 )/n 3 ] 
+ [( c 2 “ 2 c,c 3 + 2c 0 c 4 )m 4 + (c 2 - 2c 0 c 2 )m 5 + c 0 2 m 6 ]; 
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where 



m 0 = -7- {d$ tn i ~ d i n h + d \ ,n i )- 

d 7 

= { d ^ d ^ d ^ d ^) (^0^(3 ^1 ^2) (^0^7 -\- d 2 d $ ^(3^4)* 

m 2 = (^0^7 “ 44)( “^0^5 + ^1^4) + (4)^3 “ d ] d 2 ){d 2 d 1 — d 3 d 6 ) y 

m 3 ~ “ (^ 0^7 “ ^ 1 ^ 5 ) + (4)^3 “ ^1 d 2 ){d$d 2 — ^ 5 ^ 5 ), 

m 4 = — j- (^2^3 — ^4^2 + ^6 m i)» 

“0 

m s = -j- {d 2 m± - </ 4 m 3 + d 6 m 2 ) y 

a o 

m 6 = y— {d 2 ??i$ — </ 4 m 4 + d 6 m 3 ), and 

d 0 

A 7 = ^o(^l m 6 ” <^5 + ^5 m 4 “ ^7 m 3)- 

/ 8 = -77- [c 7 m 0 + (cj - 2 c 3 c 1 )m x + (c 5 2 - 2c 4 c 6 + 2c 3 c 7 )m 2 + (cl - 2 c 3 c 5 + 2 c 2 c 6 - 2c,c 7 )m 3 ] 
+ 7^- [(c 3 - 2 c 2 c 4 + 2c,c 5 - 2 c 0 c 6 )m 4 + (c 2 - 2c,c 3 + 2 c 0 c 4 )m 5 + (c 2 - 2c 0 c 2 )//z 6 + Co//j 7 ]; 

where 

m 0 = ~~T ( d 6 m l ~ d <i m 2 + d 2 m 3 ~ d o n h\ 

“8 

m \ = i d o d i - d 2 d s)( d Q d \ d 7 + d 0 d 3 d 5 + 2 d 7 d 6 ) + (d 3 d-, - dj)(dld 5 + d x d]) 

~h d^ d 3 d^ ( d 0 d 3 — d\ ^2) — d\ dgi^dfid^ d] d^) 

+ ( — d 2 d 2 + d 2 d 6 — d^d^d^d] + d\d A ) — d l d 6 (d]d 6 4 - 3 d 0 d 3 d s ) 

— c/j d 2 d 3 {d 3 d^ — d^d 3 ) + 2 d fid} d^d 3 , 

m 2 = (d 0 d 3 ~ d^d^d^d-j ~ d\ d$d% — d x d^d 2 + d 2 d 5 d 2 ) 

+ (d 3 d% — d 4 d 2 )( —d 0 d ] d 5 + d 0 d 3 — ^^2^3 + d^d^) 

— dQd s d 2 (d 0 d s — d \d^) + d ] d§(d 0 d- — d^d^), 
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2 2 2 
n h = d\ {d\d% — d 2 d 2 ) 4- ( —d 5 d$ 4- d 6 d-j)(d Q d x d 5 — d Q d 3 4- d x d 2 d 3 — d x d^) 

+ -d 0 d 5 + d x d 4 + d 2 d 3 ) - 2 d 0 d l d 1 d s , 

m 4 = ( d$dg 4- d 6 d 2 )( 2 d 0 d x d-i — d 0 d 3 d 5 4- d ] d 2 d 5 — d x d 6 ) 

+ ( — d 3 d s 4 - d A d-j)(d 0 d 3 d 2 — d x d 2 d 7 4 - ^ ^ g ) — d Q d 7 , 

m 5 = --- (^ 2 m 4 - ^ 4 a?? 3 4- d 6 m 2 — 
a Q 



m 6 = (</ 2 "* 5 - d A m A 4- rf 6 m 3 - ^ s m 2 ), 

a o 



m-j = -j- {d 2 m 6 — d A m s 4- J 6 /n 4 — d s m 3 ), and 
"o 



= 4)W'"7 - ^3 m 6 + ^5 /m 5 - d-jf72^). 



h = 3^- [ c 8 m 0 + ( c 7 - 2c 6 c 8 )m, + (c- 6 2 - 2c 5 c 7 + 2 c 4 c 8 )m 2 + (cj - 2 c 4 c 6 + 2 c 3 c 7 - 2c 2 c g )w 3 ] 
+ [(c 4 - 2 c 3 c 5 + 2c 2 c 6 - 2c,c 7 + 2c 0 c 8 )/n 4 + ( c 3 - 2 c 2 c 4 + 2c,c 5 _2c 0 c 6 )m 5 ] 

+ [(c 2 — 2 c ] c 3 + 2 c 0 c 4 )m 6 + (c 2 — 2 c 0 c 2 )m 7 + Cq/Mj]; 



where 



m Q = (d 1 m l — d 5 m 2 4- d 3 m 3 — d x m 4 ), 
“9 



W2 X — £?l(£?l£?IQ Cl 2 Cl 9 4" @ 3 &6 4" 4" 2<7 4 (2 6 £2$ d^Q-j # 7 ) 

+ £2 2 ( — ^2^6 — ^3^7 "h ^4^5 “I" 2(2 4 £? 2 ) — £? 4 , 



m 2 = aj(a 3 £29 4- a 4 £?9 — a s a s a 6 a i ~~ a i a &) & 2 { — a 2 a 9 4- tf 4 a g + < 27 ) — ^ 4 ^ 7 , 



m 3 = a](a 3 a 10 + a 4 a 10 + a 7 a 9 - a 8 ) + a 2 ( - a 2 a 10 + a 7 a 8 ) - « 4 a 7 , 



^4 — a i( a 5 a io + 2a 7 a 10 — ^ 8 ^ 9 ) + ^ 2(^7 ^9 — fl 4 a io) — °7> 



m 5 = -J- (</ 2 m 4 - d 4 m 3 + d 6 m 2 - ^m,). 

“0 
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m 6 = -j- (^ 2 m 5 - d 4 m A + d 6 m 3 - d a m 2 ), 

“o 

w 7 = -j- (d 2 m 6 - d 4 m 5 + </ 6 m 4 - </ 8 wi 3 ), 

“o 

m g = -j- (d 2 tn 2 — d A m 6 + d 6 m 5 — d & m 4 ), and 

“o 

Ag = d 0 (d ] m s — d^mj + d 5 m 6 — d 1 m 5 + </ 9 m 4 ), 

where 

a i d^d 2 d 0 d 3 

a 2 ~ d^d^ — d§d 3 

a 3 = ^3^4 “ d 2 d 5 

o.^ = d^dfr d^d-j 

a s = ^3^6 ~ d 2 d- t 

cl§ = d^d^ — d^d-j 

a i = d Y d% — d 0 dg 

a s = d 3 d% — d 2 d 9 
dg — d^dg — d$dg 

a io = ^7 d% — d 6 dg. 
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APPENDIX B. PERFORMANCE INDEX PROGRAM. 



This appendix contains a listing file of the FORTRAN source code for the per- 
formance index program. This program is used to plot the performance index curves for 
desired values. 

The program consists of main program and subroutines. The subroutine SIZE per- 
mits adjustment of plot size to both the screen and the printer. The subroutine LDFILE 
permits program input of a pre-existing data file. Subroutine MON'PRT allows selection 
of the graphics output device, be it the monitor or a wide variety of printers. 



* * * * * * ARRAY definition * * * * * * 



CHARACTER* 1 CHAR,CHG 
REAL Z( 5 1 , 5 1) 

REAL XP(202) ,YP(202) ,ZP(200) 

REAL ZPI J( 200) 

INTEGER KNXT( 200) , ITEXT( 1) 

INTEGER J,K,SEL 
REAL ZLEV( 50) 

INTEGER LDIG(50) ,LWGT(50) 

C0MM0N/PL0TTR/ I0P0RT, MODEL 

* * * * * * VARIABLE INITIAL VALUE * * * * * * 

DATA XMAX/11. 0/,YMAX/8. 0/ 

DATA XLPL0T/1. 0/,YLPL0T/0. 5/,XHPL0T/8. 0/,YHPL0T/5. 5/ 
DATA CAY/5. 0/,NRNG/3/ 

DATA HGT/. 12/ ,NDIV/4/ ,NARC/8/ 

DATA NSM/4/ 



* * * * * * ENTER INPUT PARAMETERS * * * * * * 



10 CONTINUE 
CALL CLR 
CALL CUP(1,0) 

M0DEL=99 
I0P0RT=99 
13 WRITE(* , 297 ) 

WRITE(*,305) 

READ(* ,*) NLEV 

IF (NLEV. LT. 0.0R. NLEV. GT. 50) GO TO 13 
DO 22 K=1 ,NLEV 
WRITE(*,304) K 

READ(* ,*) ZLEV(K) , LWGT(K) ,LDIG(K) 

22 CONTINUE 

WRITE(*, 300) 

READ(* ,*) NPTS 
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WRITE (*,29 7) 



* * * * * * BEGIN PROCESSING * * * * * * 



FACT=0. 9 

CALL SIZE (FACT) 

25 CONTINUE 

CALL LDFILE(NPTS ,XP , YP , ZP) 

C. .. INITIALIZE PLOT88 

CALL PLOTS(0,IOPORT, MODEL) 

C. . . DEFINE OUTPUT WINDOW 

CALL WINDOW( 0. ,0. ,XMAX,YMAX) 

C. . . REORIGIN THE ENTIRE DRAWING 
CALL PLOT(l. 0,0. 5,-3) 

XPMAX=XP( 1) 

YPMAX=YP( 1) 

XPMIN=XP( 1) 

YPMIN=YP( 1) 

DO 3 K=2 , NPTS 

IF (XP(K) . GT. XPMAX) THEN 
XPMAX=XP(K) 

END IF 

IF (YP(K). GT. YPMAX) THEN 
YPMAX=YP(K) 

END IF 

IF (XP(K).LT. XPMIN) THEN 
XPMIN=XP(K) 

END IF 

IF (YP(K).LT. YPMIN) THEN 
YPMIN=YP(K) 

END IF 

3 CONTINUE 

DO 4 J=1 ,NPTS 

XP ( J ) =XP( J ) *XHPLOT/XPMAX 
YP( J)=YP( J)*YHPLOT/ YPMAX 
XHIGH=XHPLOT 
YHIGH=YHPLOT 
XLOW=XPMIN*XHPLOT/XPMAX 
YLOW=YPM I N* YHP LOT / YPMAX 

4 CONTINUE 

C. . . ZERO GRID ARRAY Z 
DO 200 1=1,51 
DO 200 J=1 ,51 
Z(I,J)=0. 0 
200 CONTINUE 



C. . . GRID SIZE PARAMETERS 
NX=8 
NY=6 
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DX=( XPMAX-XPMIN) /FLOAT( NX- 1 ) 

DY=( YPMAX-YPMIN) /FLOAT( NY - 1 ) 

CALL FACTOR (FACT) 

CALL STAXIS ( . 13, . 20, . 15 .0. 1,2) 

CALL AXIS( 1. 0,0. 5, 'Alpha' ,-5,-7. ,0. ,XPMIN,DX) 

CALL AXIS(1. 0,0. 5, 'Beta' ,4,-5. 0,90. , YPMIN.DY) 

CALL ZGRID( Z , 51 , 51 ,NX,NY ,XLOW, YLOW,XHIGH, YHIGH,XP , YP , ZP , 

1 NPTS , CAY , NRNG , ZPI J , KNXT) 

C. . POST THE INPUT POINTS. 

CALL ZPOST( XLOW , YLOW , XHIGH , YHIGH . XLPLOT , YLPLOT , XHPLOT , YHPLOT , 

1 XP, YP, ZP ,NPTS ,2,(lx,f5. 2) ',.1389) 

IRTN=590 

C PUT A BOX AROUND THE MAP 

580 CALL PLOT( XLPLOT, YLPLOT, 3) 

CALL PLOT( XHPLOT, YLPLOT, 2) 

CALL PLOT( XHPLOT, YHPLOT, 2) 

CALL PLOT( XLPLOT, YHPLOT, 2) 

CALL PLOT( XLPLOT, YLPLOT, 2) 

IF(IRTN. EQ. 900) GOTO 900 

590 CALL PLOT(l. 0,0. 5,-999) 

CALL PL0T(1. 0,0. 5,-3) 

CALL FACTOR (FACT) 

CALL STAXIS ( . 13 , . 20 , . 15 , 0. 1 , 2) 

CALL AXIS( 1. 0,0. 5, 'Alpha' ,-5,-7. ,0. ,XPMIN,DX) 

CALL AXIS(1. 0,0.5, 'Beta' ,4, -5. 0,90. ,YPMIN,DY) 

C SMOOTH THE DATA BEFORE DRAWING 

CALL ZSMTH(Z,51,51,NX,NY,NSM) 

CALL FACTOR (FACT) 

CALL STAXIS ( . 13 , . 20 , . 15 .0. 1 , 2) 

CALL AXIS(1. 0,0. 5, 'Alpha' ,-5,-7. ,0. ,XPMIN,DX) 

CALL AXIS(1. 0,0. 5, 'Beta' ,4,-5. 0,90. , YPMIN,DY) 

CALL ZCSEG( Z , 5 1 , 5 1 , NX , NY , XLPLOT , YLPLOT , XHPLOT , YHPLOT , 

1 ZLEV , LDI G , LWGT , NLEV , HGT , ND I V , NARC ) 

IRTN=900 
GO TO 580 

900 CALL PLOT( 1. 0,1. 0,999) 

WRITE(*,*)' WOULD YOU LIKE HARDCOPY OF THESE PLOTS? (Y or N) ' 
READ(*, 401) CHAR 

IF((CHAR. EQ. ' Y' ). OR. (CHAR. EQ. 'y')) THEN 
CALL MONPRT 
GOTO 25 
ENDIF 

WRITE(*,*)’ WOULD YOU LIKE TO DO ANOTHER ? (Y or N)' 
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READ(*,401) CHAR 

IF(( CHAR. EQ. ' Y * ). OR. (CHAR. EQ. 'y' )) GOTO 10 



297 FORMAT(lx) 



300 FORMATC ' 

304 FORMATC' 

305 FORMATC ' 
401 FORMATC A) 



3. Enter number of data file point: 
2. Enter PI values , lwgt,ldig # ', 
l.How many constant PI curves # = 




C. . . end of program 



STOP 

END 



C 

C 

C 

C 

C 



SUBROUTINE CLR 

SUBROUTINE TO ERASE THE DISPLAY. 
CALL , CALL CLR 



* 

* 



SUBROUTINE CLR 
CHARACTER* 1 C1,C2,C3,C4 
INTEGER*2 ICC 4) 

EQUIVALENCE C Cl , IC( 1) ) , ( C2, IC( 2) ) , C C3 , ICC 3) ) , CC4 , IC(4) ) 
DATA IC/16#1B,16//5B, 16#32 ,16#4A/ 

C. . . WRITE ESCAPE CODE TO DISPLAY 
WRITE(* , 1) C1,C2,C3,C4 
1 FORMATC lx, 4al) 

RETURN 

END 



C MOVE THE CURSOR TO LINE N, COLUMN M. 
C CALL, CALL CUP(N,M) 

C N - NEW CURSOR LINE 

C M - NEW CURSOR COLUMN 

C - 



SUBROUTINE CUP CN,M) 

CHARACTER* 1 C1,C2,C5,C8,LC(5) 

CHARACTER*5 CBUFF 
INTEGER* 2 IC(4) 

EQUIVALENCE ( Cl , IC( 1) ) , ( C2 , ICC 2) ) ,CC5 , IC( 3) ) , (C8 , IC(4) ) , 
*( CBUFF, LC( 1)) 

DATA IC/16#1B , 16#5B , 16#3B , 16#66/ 

L=10000+100*N+M 

WRITE(CBUFF,2)L 
2 FORMATC i5) 

WRITEC*, 1) C1,C2,LCC2),LC(3),C5,LCC4),LC(5),C8 
1 FORMATC lx, 8a 1, ) 

RETURN 



C 

C 



SUBROUTINE CUP 



END 



C 

C 
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nnnnnroh^h-* non 



SUBROUTINE LDFILE -- ALLOWS LOADING OF DATA FROM AN EXISTING FILE * 

“V 



SUBROUTINE LDFILE (NPTS ,XP , YP , ZP) 

REAL XP( 202) , YP( 202) , ZP( 200) 

INTEGER NPTS, I 
CHARACTER* 12 FILEIN,ANS*1 
LOGICAL EXIST 
1 CALL CLR 
WRITE(*,100) 

100 F0RMAT( IX, //// . 10X, ' 

* ' ,/.10X, *** THIS IS THE ROUTINE TO LOAD DATA FROM A FILE ** 

**' ,/,10X, f ' 

*,//,10X, 'FILE NAME SHOULD NOT EXCEED 8 CHARACTERS' / , 

*10X, ' FILE EXTENSION SHOULD NOT EXCEED 3 CHARACTERS*) 

WRITE(*, 101) 

101 FORMAT( IX, //// , 10X, 'WHAT IS THE FILE NAME (FN) AND EXTENSION (EXT) 
*? ' ,//,10X,' ENTER IN FORM FN. EXT (E.G. PARAPLAN. INP) ==> ') 

CALL CUP( 18,57) 

READ(*,102) FILEIN 

102 FORMAT( A12) 

I NQU I RE ( F I LE=F I LE I N , EXI ST=EX I ST ) 

IF (EXIST) THEN 

OPEN( 7 , FILE=FILEIN , STATUS= ' OLD ' , ACCESS=' SEQUENTIAL ' ) 

v T cr 

WRITE(* , 103) 

103 F0RMAT(///,10X, 'THAT FILE DOES NOT EXIST/, 10X, 

* 'DO YOU WANT TO TRY AGAIN? (Y/N)') 

READ(*, ' ( Al) ' ) ANS 

IF( ( ANS .EQ. ’ Y’ ) .OR. (ANS . EQ. 'y')) THEN 
GOTO 1 
ELSE 
GOTO 2 
END IF 
END IF 

DO 190 I = l.NPTS 

READ( 7 , 120) XP( I ) , YP( I ) , ZP( I ) 

20 FORMAT( 1X,F10. 2 , 3X,F10. 2 , 3X,F10. 3) 

90 CONTINUE 

CLOSE( 7 , STATUS=' KEEP 1 ) 

RETURN 

END 



* 

SUBROUTINE MONPRT -- ALLOWS DISPLAY ON MONITOR OR PRINTER * 

* 



SUBROUTINE MONPRT 
INTEGER IOPORT, MODEL, SEL 
COMMON/ PLOTTR/ IOPORT, MODEL 

CALL CLR 



62 



80 

81 

87 

82 

83 



WRITE(* , 80) 

WRITE(* ,81) 

WRITE(* , 82) 
WRITE(*,81) 

WRITE (*, 83) 

WRITE(*, 87) 

WRITE (*,84) 
WRITE(*,81) 
WRITE(*,93) 
FORMAT(lx,///,10x f '■ 



, / , lOx, ' ' ,21x, 'PRINTER/OUTPUT MENU' ,22x, ' ') 



*. 

FORMAT( 10x, 1 

* ') 

FORMAT( lOx , ' ++++++++++++ ++F H I I II I II H H I + 1 I H + f + R +-I-++++++++- 
*++++++++++ ' ) 



FORMAT( lOx, 1 
FORMAT( lOx, ' 
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,2x, 'PRINTER NO. ' ,2x, ' ' ,20x, ' PRINTER' , 19x, ' ') 
, 3x, ' 1 ' , 2x, ' ' ,4x, 'Epson FX-80, All' 



* > / > lOx, ' ' 
* , 25x , ' 


,3x, ' 


2 


' , 2x, ' ' 


^x/Epson FX-100, All 1 


*> / » lOx, ' ' 
*»25x, ' ' 


,3x, ' 


3 


t 0 It 

,2x. 


,4x, ! Epson MX-100, All 1 


* , / , 1 Ox , ' 
* , 26x , ' ' 


X 

CO 


4 


' , 2x, ' ' 


,4x, ! Epson RX-80, All 1 


*,/ » lOx, ' ’ 


, 3x, 


5 


' , 2x, ' ' 


^x^Epson MX-80 6c IBM Printer 1 


*,/, lOx, ' ' 
* , 1 7 x , ' ' 


> 3x, ' 


6 


' , 2x, ' ' 


,4x, ! HP 7470A Graphics Plotter 1 


* , / , 1 Ox , ' ' 
* > 1 7 x , ' ' 


, 3x, ' 


7 


' ,2x, ' ' 


,4x, f HP 7475A Graphics Plotter 1 


* , / , 1 Ox , ' ' 
*> 18x, ' ' 


,3x, ' 


8 


' , 2x, ’ ' 


,4x, ! HP 758xB Series Plotters 1 


* , / , lOx , ' ’ 
* , 24x , ' ') 
format( lOx 


,3x, ' 

,' ' ,3x, 


9 


' , 2x , ' ' 


,4x, ! HP 2686A Laser Jet 1 


*' 10 
*» I6x, ' ' 


' , 2x, ' 


' ,4x, 


' Graphics 


Monitor (default) 1 


*,/> lOx, ' ' 
* , 19x , ' ' 


,3x, ' 


11 


' , 2x , ’ ' 


, 4x , 1 HARDWARE Interface Menu 1 


*>/ > lOx, ' ,3x, ' 

*RT and MODEL' 

* , 2x , ' ' 


12 


’ 2x ’ ' 


,4x, f Input PLOT88 Values for IOPO 


* , / , 1 Ox , ' 
^ 


>3x. 


99 


1 o ft 

,2x, 


, 4x , 1 RETURN to LOAD Data’,23x, 
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FORMAT( lx, ///,15x, 'ENTER INTEGER NUMBER FOR SELECTION ==> ') 
CALL CUP( 24,56) 

READ(*,*) SEL 
IF (SEL .EQ. 1) THEN 
IOPORT = 0 
MODEL = 5 

ELSE IF (SEL . EQ. 2) THEN 
IOPORT = 0 
MODEL = 15 

ELSE IF (SEL . EQ. 3) THEN 
IOPORT = 0 
MODEL = 11 
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o o o o o 



ELSE IF (SEL . EQ. 
IOPORT = 0 
MODEL = 1 


4) 


THEN 


ELSE IF (SEL . EQ. 
IOPORT = 0 
MODEL = 1 


5) 


THEN 


ELSE IF (SEL . EQ. 
IOPORT = 9600 
MODEL =20 


6) 


THEN 


ELSE IF (SEL . EQ. 
IOPORT = 9600 
MODEL = 30 


7) 


THEN 


ELSE IF (SEL . EQ. 
IOPORT = 9600 
MODEL = 80 


8) 


THEN 


ELSE IF (SEL . EQ. 
IOPORT = 9600 
MODEL = 60 


9) 


THEN 


ELSE IF (SEL . EQ. 
CALL PORT 


11) 


THEN 


ELSE IF (SEL . EQ. 


12) 


THEN 



CALL CLR 
WRITE(*, 100) 

100 F0RMAT( lx,//// ,5x, ' INPUT VALUE TO BE USED BY THE PL0T88 GRAPHICS 
* PACKAGE* ,/,5x, 'FOR IOPORT =>',//) 

CALL CUP( 7,21) 

READ(*,*) IOPORT 
WRITE (*,110) 

110 FORMAT( lx, ///, 5x, ' INPUT VALUE TO BE USED BY THE PL0T88 GRAPHICS 
♦PACKAGE' ,/,5x,' FOR MODEL =>',//) 

CALL CUP( 12,20) 

RE AD(*,*) MODEL 
ELSE 

IOPORT = 99 
MODEL =99 
END IF 

WRITE(*,98) IOPORT, MODEL 

98 FORMAT( / / , lOx , ' PL0T88 WILL USE THESE VALUES TO OUTPUT GRAPHICS: ' 
*,//,10x, ' IOPORT = ' ,14, lOx, 'MODEL = ’,14,///) 

PAUSE 

RETURN 

END 



* 

SUBROUTINE PORT -- ALLOWS SELECTION OF OUTPUT PRINTER PORT * 

* 



SUBROUTINE PORT 

INTEGER IOPORT, MODEL, SEL 

COMMON PLOTTR/ IOPORT, MODEL 

CALL CLR 
WRITE(*, 79) 

79 FORMATC lx,///, 15x,' SELECTION OF A PRINTER IN THE PREVIOUS MENU AUT 
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*OMATICALLY’ ,/,10x, 'SELECTS THE MOST COMMONLY ASSOCIATED OUTPUT POR 
*T, BE IT PARALLEL ',/,10x, 

*' (LPT) OR SERIAL, FOR THAT PARTICULAR DEVICE. IF THE PROGRAM GRAPH 
*ICS' , / , lOx, 

* ' ARE NOT BEING OUTPUT CORRECTLY, USE THIS MENU TO PROPERLY ROUTE T 
*HE ' ,/ , lOx, 

*' GRAPHICS DATA TO THE OUTPUT DEVICE. ’,////) 

PAUSE 



CALL CLR 
WRITE(*,80) 

WRITE(*,81) 

WRITE(*,82) 

WRITE(*,81) 

WRITE(*,83) 

WRITE(*,87) 

WRITE (*,84) 

WRITE(*,81) 

WRITE(*,93) 

80 FORMAT(lx,///,10x,’ 

* 1 ,/ , lOx, ' ’ ,19x, 'HARDWARE INTERFACE MENU' 

*, 20 x, ! ! ) 

81 FORMAT( lOx, ' 

* ') 

8 7 FORMAT( lOx , ' ^++++++^+++++-H-++++++++++4 +++++-H-++-H--H-+-H-+++-H--H-+++ 

*++++++++++' ) 



82 FORMAT( lOx, ’, 2x, ’ SELECT NO. ', 3x , ' , 20x , PORT , 22x, ’) 

83 FORMAT( lOx, ' ’,3x,’ 1 ',2x,' 14x, ' LPT1 Printer Port' 





*,15x,’ ' 
* , / , lOx, ' 


» 3x, ' 


2 


' , 2 x, ' 


' ,14x, 'LPT2 


Printer Port' 




*, 15x, ' 

*> / > lOx, ' ' 


,3x, ' 


3 


' > 2 x, ' 


' , 14x, ' LPT3 


Printer Port' 




*,15x,' ' 
*>/ , lOx, ’ 


>3x, ' 


4 


' > 2 x, ' 


' , 14x, ' C0M1 


Serial Port' 




*,16x,’ ' 

*,/, 10 x,' ’ 


,3x,' 


5 


' > 2 x, ' 


' ,14x, 'COM2 


Serial Port' 


84 


* I 6 x ' ' 1 

FORMAT( lOx 
*' 99 


,' ' ,3x, 
' > 2 x , ' 


' , 14x, ' 


1 RETURN 


to LOAD Data' ,13x,' ’) 


93 


FORMAT( lx , 


III , 15x, 


' ENTER 


INTEGER 


NUMBER FOR 


SELECTION => 




CALL CUP( 21,56) 
READ(*,*) SEL 
IF (SEL .EQ. 1) THEN 
IOPORT = 0 









ELSE IF (SEL . EQ. 2) THEN 
IOPORT = 2 

ELSE IF (SEL . EQ. 3) THEN 
IOPORT = 3 

ELSE IF (SEL . EQ. 4) THEN 
IOPORT = 0 

ELSE IF (SEL . EQ. 5) THEN 
IOPORT =50 
ELSE 

IOPORT =99 
ENDIF 

IF ((SEL .EQ. 4) .OR. (SEL . EQ. 5)) THEN 
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60 



61 

62 

63 



CALL CLR 
WRITE(*,60) 

WRITE (*,61) 

WRITE(*,62) 

WRITE(*,61) 

WRITE (*, 63) 

WRITE(*,61) 

WRITE(*,93) 

FORMATC lx, /// , lOx . ' 

* , / , lOx, ' ' , 16x, ' BAUD (data transfer) RATE MENU' 

*> 16x, ' ') 

FORMAT(10x ' 

* ' ) 

FORMAT( lOx, ' ' ,2x, 'SELECT NO. ’ ,3x,' ', 15x, ' BAUD RATE ', 22x, ' ') 

FORMAT( lOx, ' ' , 3x, ' 1 ',2x,' ',18x,' 300' 

* , 24x , ' ' 



*,/ , lOx, ' 
*, 24x, ' 


' , 3x, ' 


2 


X 

CnI 


' ,18x, 


' 1200 


*,/ ,10x, ' 
*, 24x, ’ ' 


' > 3x , ' 


3 


' , 2x, ' 


',18x, 


'4800 


*> / , lOx, ' 


' , 3x, ' 


4 


' ,2x,' 


' ,18x. 


'9600 



* , 24x, ' ' ) 

CALl’cUP( 18,56) 

READ(* ,*) SEL 
IF (SEL .EQ. 1) THEN 
IOPORT = IOPORT + 300 
ELSE IF (SEL . EQ. 2) THEN 
IOPORT = IOPORT + 1200 
ELSE IF (SEL . EQ. 3) THEN 
IOPORT = IOPORT + 4800 
ELSE 

IOPORT = IOPORT + 9600 
ENDIF 



50 



51 

52 

53 



CALL CLR 
WRITE(*,50) 

WRITE(* ,5 1) 

WRITE(*,52) 

WRITE (*,51) 

WRITE(* ,53) 

WRITE(*,51) 

WRITE(*,93) 

FORMAT( lx , /// , lOx , ' 

* ' ,/ , lOx, ' ' ,25x, 'PARITY MENU' 

* , 26x, ' ’) 

FORMAT( lOx, ' 

* ') 



FORMAT( lOx, ' ' ,2x, 'SELECT NO. ' ,3x, ' ', 17x, ' PARITY' ,23x, 
FORMAT( lOx, ' ' , 3x, ' 1 ',2x,' ',14x,' NO Parity' 



*,22x, ' 
*, / , lOx, 



’ >3x, ' 



' , 2x , ' ' ,14x, 'EVEN Parity' 



* j 2 lx , 
*>/ j lOx, 
* , 22x, ' 



) 



' , 3x, ' 



CALL CUP( 17,56) 



' , 2x, ' ' ,14x, 'ODD Parity' 



’) 
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o o o n o 



READ(*,*) SEL 
END IF 

IF (SEL .EQ. 3) THEN 
IOPORT = IOPORT + 1 
ELSE IF (SEL . EQ. 2) THEN 
IOPORT = IOPORT + 2 
EENDIF 
99 RETURN 
END 



* 

SUBROUTINE SIZE -- ALLOWS ADJUST OF OUTPUT PLOT * 

* 



SUBROUTINE SIZE (FACT) 

REAL FACT 

C. . . ADJUST SIZE OF OUTPUT PLOT 
CALL CLR 
WRITE(*, 15) 

15 FORMAT( lx, //// , lOx , 1 

* 1 ,/,10x,'*** This is the routine to Adjust the plot size ** 

**',/, 10x, f ' 

*111, lOx, ' Output plots are currently sized to fill the console sere 

*en. , 

* lOx, Would you like to adjust the plot size? ,/,10x, 

- Enter y or n => ) 

CALL CUP( 14,31) 

READ(* , 1 ( A) ' ) CHG 

IF ((CHG .EQ. 'n') .OR. (CHG . EQ. 'N')) THEN 
GO TO 17 
END IF 
WRITE(*, 16) 

16 FORMAT( III , lOx , ' A default plot factor size of 0.9 is used to fill 
*the console screen. lOx, 

*'This plot size will be halved by entering a value of 0.45. 

*10x , ' It will be doubled with a factor entry of 1.8. ',//,10x, 

*' Enter factor size ==> ') 

CALL CUP( 22,32) 

READ(*,*) FACT 

17 RETURN 
END 
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APPENDIX C. D.S.L SIMULATION PROGRAM AND DATA FILE. 



This appendix contains three listing files and its resulting data files for the D.S.L 
simulation programs used in examples 4-1, 4-2 and 4-3. These programs are used to 
produce the data files which generate the plots of the performance index curves on the 
parameter plane. 



A. SIMULATION PROGRAM FOR EXAMPLE 4-1. 

TITLE PERFORMANCE INDEX CALCULATION 
ARRAY Z(l) ,P(4) ,C(2) ,D(2) 

TABLE Z(l) = 100. ,P(1-4) = 1. 0,11. 0,10. 0,0. 0,D(1)=1. 0 

PARAM A=0. 0,B=0. 0,Y0=0. 0 

INITIAL 

C(1)=A 
C(2)=B 
D( 2)=B 
BETA=B 
ALPHA=A 
DYNAMIC 

R=STEP( 0. ) 

DERIVATIVE 

Y1=TRNFR( 0,3,Y0,Z,P,E) 

M=TRNFR( 1,1,Y0,C,D,Y1) 

E=R-M 

ABE=ABS(E) 

PI1=INTGRL( Y0 ,E*E) 

PI2=INTGRL( Y0 , ABE*TIME ) 

PI3=INTGRL( YO , ABE) 

CONTROL FINTIM=15 ,DELT=0. 01 
PRINT 15. , ALPHA, BETA, PI1,PI2,PI3 
END 

PARAM A=4. ,B=10. 

END 

PARAM A=5. , B=7. 5 
END 

PARAM A=5. ,B=10. 

END 

PARAM A=5. ,B=15 
END 

PARAM A=6. ,B=10 
END 

PARAM A=7. , B=30 
END 

PARAM A=7. , B=40 
END 

PARAM A=7. 5 ,B=15 
END 

PARAM A=7. 5 , B=20. 
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END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 



A=7. 5 ,B=25. 
A=7. 5 ,B=35 
A=7. 5 ,B=45. 
A=7. 5,B=50. 
A=10. ,B=10. 
A=10. ,B=15 
A=10. ,B=20. 
A=10. , B=25. 
A=10. ,B=35 
A=10. ,B=40. 
A=12. 5 , B=10. 
A=12. 5 ,B=15 
A=12. 5 ,B=20. 
A=12. 5 ,B=25. 
A=12. 5 , B=35 
A=12. 5 , B=45. 
A=12. 5 ,B=50. 
A=15. ,B=10. 
A=15. ,B=15 
A=15. ,B=20. 
A=15. ,B=25. 
A=15. ,B=30. 
A=15. ,B=35 
A=15. , B=40. 
A=15. ,B=45 
A=15. ,B=50. 
A=17. 5 ,B=15 
A=17. 5 ,B=20. 



END 

PARAM A=17.5,B=25 
END 

PARAM A=17. 5 ,B=35 
END 

PARAM A=17. 5 , B=40 
END 

PARAM A=20. ,B=15 
END 

PARAM A=20. ,B=20. 
END 

PARAM A=20. ,B=25. 
END 

PARAM A=20. ,B=35 
END 

PARAM A=20. , B=40. 
END 

PARAM A=20. ,B=30. 
END 

PARAM A=20. ,B=60. 
END 

PARAM A=20. ,B=90 
END 

PARAM A=25. ,B=15 
END 

PARAM A=25. ,B=20. 
END 

PARAM A=25. ,B=25. 
END 

PARAM A=25. ,B=35 
END 

PARAM A=25. ,B=40. 
END 

PARAM A=25. ,B=30. 
END 

PARAM A=30. ,B=20. 
END 

PARAM A=30. ,B=25. 
END 

PARAM A=30. , B=30 
END 

PARAM A=30. ,B=35. 
END 

PARAM A=30. , B=40. 
END 

PARAM A=50. , B=20. 
END 

PARAM A=50. ,B=25. 
END 

PARAM A=50. , B=30 
END 

PARAM A=50. ,B=35. 
END 

PARAM A=50. , B=40. 
END 

PARAM A=50. ,B=50. 
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END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

STOP 



A=50. , B 
A=60. ,B 
A=60. ,B 
A=60. ,B 
A=80. ,B 
A=80. ,B 
A=80. ,B 
A=80. , B 
A=100. , 
A=100. , 
A=100. , 
A=100. , 
A=100. , 
A=100. , 
A=120. , 
A=120. , 
A=130. , 
A=130. , 
A=140. 
A=140. 



= 60 . 

=30. 

=70 

=90. 

=30. 

=50. 

=70. 

=90. 

B=30. 

B=40. 

B=50. 

B=60. 

B=80. 

B=100 

B=35. 

B=50. 

B=35. 

B=50. 

,B=35. 

, B=50. 
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B. DATA FILE FOR EXAMPLE 4-1. 



ALPHA 

0.0000 

4.0000 

5.0000 

5.0000 

5.0000 

6.0000 

7.0000 

7.0000 

7.5000 

7.5000 

7.5000 

7.5000 

7.5000 

7.5000 

10.000 

10.000 

10.000 

10.000 

10.000 

10.000 

12.500 

12.500 

12.500 

12.500 

12.500 

12.500 

15.000 

15.000 

15.000 

15.000 

15.000 



BETA 

0.000 

10.000 

7.500 

10.000 

15.000 

10.000 

30.000 

40.000 

15.000 

20.000 

25.000 

35.000 

45.000 

50.000 

10.000 

15.000 

20.000 

25.000 

35.000 

40.000 

10.000 

15.000 

20.000 

25.000 

45.000 

50.000 

10.000 

15.000 

20.000 

25.000 

30.000 



IES 

1.00000 

0.32856 

0.30910 

0.2S481 

0.30855 

0.26177 

0.33793 

0.413S7 

0.23340 

0.25469 

0.2S477 

0.35289 

0.42382 

0.45935 

0.25001 

0.20401 

0.21075 

0.22905 

0.27603 

0.30143 

0.28031 

0.19186 

0.18744 

0.19S06 

0.27154 

0.29198 

0.34726 

0.18889 

0.17409 

0.17SS5 

0.1S996 



IAET 

1.50000 

0.47462 

0.53308 

0.33893 

0.42223 

0.30259 

0.66355 

1.1865 

0.1S931 

0.25123 

0.3S517 

0.76717 

1.27460 

1.57120 

0.42101 

0.15S25 

0.14365 

0.19092 

0.3S136 

0.505S7 

0.6667S 

0.16546 

0.11826 

0.12439 

0.381S0 

0.47S37 

1.23910 

0.1S998 

0.11261 

0.10424 

0.11504 



IAE 

1.00000 

0.64140 

0.64346 

0.54590 

0.60431 

0.50S02 

0.71529 

0.92510 

0.42445 

0.47959 

0.56794 

0.7602S 

0.95456 

1.050S0 

0.54982 

0.37735 

0.37559 

0.4235S 

0.55835 

0.62925 

0.66477 

0.37150 

0.33404 

0.35079 

0.55436 

0.61059 

0.87SS9 

0.38491 

0.31716 

0.31487 

0.33643 
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15.000 

15.000 

15.000 

15.000 
17.500 
17.500 
1 7.500 
17.500 
17.500 

20.000 
20.000 
20.000 
20.000 
20.000 
20.000 
20.000 
20.000 
25.000 
25.000 
25.000 
25.000 
25.000 

25.000 

30.000 
30.000 
30.000 
30.000 

30.000 

50.000 
50.000 
50.000 
50.000 
50.000 
50.000 
50.000 



35.000 

40.000 

45.000 

50.000 

15.000 

20.000 

25.000 

35.000 

40.000 

15.000 

20.000 

25.000 

35.000 

40.000 

30.000 

60.000 

90.000 

15.000 

20.000 

25.000 

35.000 

40.000 

30.000 

20.000 

25.000 

30.000 

35.000 

40.000 

20.000 

25.000 

30.000 

35.000 

40.000 

50.000 

60.000 



0.203 S 5 

0.21918 

0.23531 

0.25193 

0.19226 

0.16644 

0.16621 

0.18434 

0.19669 

0.20112 

0.16250 

0.15766 

0.17018 

0.18022 

0.16203 

0.22763 

0.30424 

0.23759 

0.16216 

0.14800 

0.15137 

0.15796 

0.14718 

0.16949 

0.14448 

0.13 SS 9 

0.13990 

0.14392 

0.31020 

0.16519 

0.13466 

0.12415 

0.12067 

0.12200 

0.12792 



0.14462 

0.19111 

0.24730 

0.31143 

0.23148 

0.11377 

0.09545 

0.10984 

0.13387 

0.29311 

0.12059 

0.09084 

0.09490 

0.10688 

0.08883 

0.23630 

0.56953 

0.51682 

0.15589 

0.09050 

0.07941 

0.08447 

0.07807 

0.21232 

0.10924 

0.07458 

0.06935 

0.07283 

1.10170 

0.26752 

0.15155 

0.10402 

0.07715 

0.05161 

0.05862 



0.37180 
0.41479 
0.45996 
0.50594 
0.41186 
0.31237 
0.29500 
0.32730 
0.35790 
0.45121 
0.31537 
0.28397 
0.29932 
0.32112 
0.28474 
0.44709 
0.653 S 1 
0.57663 
0.33829 
0.27819 
0.26559 
0.27703 
0.26276 
0.37707 
0.28776 
0.25551 
0.24805 
0.25209 
0.81409 
0.39313 
0.29655 
0.25509 
0.2331 1 
0.21576 
0.22270 
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60.000 

60.000 

60.000 

S0.000 

50.000 

80.000 

80.000 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

120.00 

120.00 

130.00 

130.00 

140.00 

140.00 



30.000 

70.000 

90.000 

30.000 

50.000 

70.000 

90.000 

30.000 

40.000 

50.000 

60.000 

50.000 

100.000 

35.000 

50.000 

35.000 

50.000 

35.000 

50.000 



0.14256 

0.12363 

0.137S4 

0.17S46 

0.10841 

0.1093S 

0.11S09 

0.27014 

0.12869 

0.10764 

0.10224 

0.10374 

0.11050 

0.20936 

0.11045 

0.251S6 

0.11297 

0.32030 

0.11619 



0.20934 
0.054 SS 
0.07347 
0.38S18 
0.10385 
0.05337 
0.05032 
0.88639 
0.23817 
0.14872 
0.10517 
0.06032 
0.04206 
0.5S01S 
0.19635 
0.S085S 
0.22142 
1.26700 
0.24697 



0.33390 

0.21479 

0.24487 

0.45147 

0.22883 

0.19932 

0.20455 

0.71319 

0.31758 

0.24754 

0.21883 

0.19541 

0.19150 

0.54931 

0.26939 

0.66742 

0.28185 

0.S5656 

0.29524 
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C. SIMULATION PROGRAM FOR EXAMPLE 4- 



TITLE PERFORMANCE INDEX CALCULATION 
ARRAY A(1),B(4) 

TABLE A(l)=50. ,B(1)=1. 0,B(4)=0. 0 
PARAM K1=0. 2,K2=0. 0,Y0=0. 0 
INITIAL 

B(2)=3+50*K2 

B(3)=2+50*K1 

BETA=K2 

ALPHA=K1 

DYNAMIC 

R=STEP( 0. ) 

DERIVATIVE 

M=TRNFR(0,3,Y0,A,B,E) 

E=R-M 

ABE=ABS(E) 

PI 1=INTGRL( Y0 , E*E ) 

PI2=INTGRL( Y0 , ABE*TIME ) 
PI3=INTGRL( Y0 , ABE) 

CONTROL FINTIM=15 ,DELT=0. 01 
PRINT 15. , ALPHA .BETA, PI 1,PI2,PI3 
END 

PARAM K1=0. 2 ,K2=0. 1 
END 

PARAM K1=0. 2 ,K2=0. 2 
END 

PARAM K1=0. 2,K2=0. 3 
END 

PARAM K1=0. 2 ,K2=0. 4 
END 

PARAM K1=0. 2,K2=0. 5 
END 

PARAM K1=0. 2 ,K2=0. 6 
END 

PARAM K1=0. 2 ,K2=0. 7 
END 

PARAM K1=0. 2,K2=0. 8 
END 

PARAM K1=0. 2,K2=0. 9 
END 

PARAM K1=0. 2 ,K2=1. 0 
END 

PARAM K1=0. 3,K2=0. 1 
END 

PARAM K1=0. 3,K2=0. 2 
END 

PARAM K1=0. 3 ,K2=0. 3 
END 

PARAM K1=0. 3 ,K2=0. 4 
END 

PARAM K1=0. 3,K2=0. 5 
END 

PARAM K1=0. 3,K2=0. 6 
END 

PARAM K1=0. 3 ,K2=0. 7 
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END 

PARAM K1=0. 3 ,K2=0. 8 
END 

PARAM K1=0. 3 ,K2=0. 9 
END 

PARAM K1=0. 3,K2=1. 0 
END 

PARAM K1=0. 4 ,K2=0. 1 
END 

PARAM K1=0. 4 ,K2=0. 2 
END 

PARAM K1=0. 4 ,K2=0. 3 
END 

PARAM K1=0. 4 ,K2=0. 4 
END 

PARAM K1=0. 4 ,K2=0. 5 
END 

PARAM K1=0. 4,K2=0. 6 
END 

PARAM K1=0. 4 ,K2=0. 7 
END 

PARAM K1=0. 4,K2=0. 8 
END 

PARAM K1=0. 4 ,K2=0. 9 
END 

PARAM K1=0. 5 ,K2=0. 1 
END 

PARAM K1=0. 5 ,K2=0. 2 
END 

PARAM K1=0. 5 ,K2=0. 3 
END 

PARAM K1=0. 5 ,K2=0. 4 
END 

PARAM K1=0. 5 ,K2=0. 5 
END 

PARAM K1=0. 5 ,K2=0. 6 
END 

PARAM K1=0. 5 ,K2=0. 7 
END 

PARAM K1=0. 5 ,K2=0. 8 
END 

PARAM K1=0. 5 ,K2=0. 9 
END 

PARAM K1=0. 5,K2=1. 0 
END 

PARAM K1=0. 6 ,K2=0. 1 
END 

PARAM K1=0. 6 ,K2=0. 2 
END 

PARAM K1=0. 6 ,K2=0. 3 
END 

PARAM K1=0. 6 ,K2=0. A 
END 

PARAM K1=0. 6 ,K2=0. 5 
END 

PARAM K1=0. 6 ,K2=0. 6 
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END 

PARAM K1=0. 6,K2=0. 7 
END 

PARAM K1=0. 6 ,K2=0. 8 
END 

PARAM K1=0. 6 ,K2=0. 9 
END 

PARAM K1=0. 6 ,K2=1. 0 
END 

PARAM K1=0. 7 ,K2=0. 1 
END 

PARAM K1=0. 7 ,K2=0. 2 
END 

PARAM K1=0. 7 ,K2=0. 3 
END 

PARAM K1=0. 7 ,K2=0. 4 
END 

PARAM K1=0. 7 ,K2=0. 5 
END 

PARAM K1=0. 7 ,K2=0. 6 
END 

PARAM K1=0. 7 ,K2=0. 7 
END 

PARAM K1=0. 7 ,K2=0. 8 
END 

PARAM K1=0. 7 ,K2=0. 9 
END 

PARAM K1=0. 7 ,K2=1. 0 
END 

PARAM K1=0. 8 ,K2=0. 1 
END 

PARAM K1=0. 8 ,K2=0. 2 
END 

PARAM K1=0. 8 ,K2=0. 3 
END 

PARAM K1=0. 8 ,K2=0. 4 
END 

PARAM K1=0. 8 ,K2=0. 5 
END 

PARAM K1=0. 8 ,K2=0. 6 
END 

PARAM K1=0. 8 ,K2=0. 7 
END 

PARAM K1=0. 8 ,K2=0. 8 
END 

PARAM K1=0. 8 ,K2=0. 9 
END 

PARAM K1=0. 8 ,K2=1. 0 
END 

PARAM K1=0. 9 ,K2=0. 1 
END 

PARAM K1=0. 9 ,K2=0. 2 
END 

PARAM K1=0. 9 ,K2=0. 3 
END 

PARAM K1=0. 9 ,K2=0. 4 
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END 

PARAM K1=0. 9 ,K2=0. 5 
END 

PARAM K1=0. 9 ,K2=0. 6 
END 

PARAM K1=0. 9 ,K2=0. 7 
END 

PARAM K1=0. 9,K2=0. 8 
END 

PARAM K1=0. 9 ,K2=0. 9 
END 

PARAM K1=0. 9 ,K2=1. 0 
END 

PARAM Kl=l. 0 ,K2=0. 1 
END 

PARAM Kl=l. 0 ,K2=0. 2 
END 

PARAM Kl=l. 0 ,K2=0. 3 
END 

PARAM Kl=l. 0 ,K2=0. 3 
END 

PARAM Kl=l. 0 ,K2=0. 4 
END 

PARAM Kl=l. 0 ,K2=0. 5 
END 

PARAM Kl=l. 0,K2=0. 6 
END 

PARAM Kl=l. 0 ,K2=0. 7 
END 

PARAM Kl=l. 0 ,K2=0. 8 
END 

PARAM Kl=l. 0 ,K2=0. 9 
END 

PARAM Kl=l. 6 ,K2=1. 0 
END 

STOP 
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D. DATA FILE FOR EXAMPLE 4-2. 



ALPHA 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.50000 



BETA 

0.00000 

0.10000 

0.20000 

0.30000 

0.40000 

0.50000 

0.60000 

0.70000 

0.S0000 

0.90000 

1.00000 

0.10000 

0.20000 

0.30000 

0.40000 

0.50000 

0.60000 

0.70000 

O.SOOOO 

0.90000 

0.10000 

0.20000 

0.30000 

0.40000 

0.50000 

0.60000 

0.70000 

0.80000 

0.90000 

1.00000 

0.10000 



IES 

1S23.40 

0.81565 

0.91714 

1.09550 

1.28S50 

1.4S5S0 

1.67S90 

1.87450 

2.06100 

2.23130 

2.40750 

0.54210 

0.66416 

0.S02S1 

0.94561 

1.09000 

1.23450 

1.37940 

1.52230 

1.66200 

0.47397 

0.57S05 

0.6SS21 

0.S0004 

0.91257 

1.02540 

1.13840 

1.25100 

1.36320 

1.47550 

0.46277 



1AET 

1177.4 

4.9245 

6.0067 

8.5450 

11.4640 

14.6030 

17.1180 

20.2920 

22.7S20 

24.4710 

26.7370 

1.3228 

2.3429 

3.SSS4 

5.7563 

7.S763 

9.9336 

12.3130 

14.3720 

16.1900 

0.5706 

I. 2171 
2.1595 
3.3445 
4.7404 
6.2509 
7.94S4 
9.5631 

II. 1670 
13.00S0 
0.3224 



IAE 

101.49 

1.8653 

2.0S13 

2.5026 

2.9349 

3.3547 

3.3763 

4.7094 

4.3853 

4.622S 

4.8S62 

1.0492 

1.3567 

1.7132 

2.0709 

2.4222 

2.7476 

3.0747 

3.3626 

3.6201 

0.7763 

1.0504 

1.3389 

1.6271 

1.9125 

2.1895 

2.4640 

2.7191 

2.9616 

3.2087 

0.6638 
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0.50000 

0.50000 

0.50000 

0.50000 

0.50000 

0.50000 

0.50000 

0.50000 

0.50000 

0.60000 

0.60000 

0.60000 

0.60000 

0.60000 

0.60000 

0.60000 

0.60000 

0.60000 

0.60000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.S0000 

o.soooo 

0.80000 

O.SOOOO 

o.soooo 

o.soooo 



0.20000 


0.55073 


0.30000 


0.64156 


0.40000 


0.73322 


0.50000 


0.82524 


0.60000 


0.91744 


0.70000 


1.00970 


o.soooo 


1.10200 


0.90000 


1.19420 


1.00000 


1.2S650 


0.10000 


0.47534 


0.20000 


0.55087 


0.30000 


0.6279S 


0.40000 


0.70557 


0.50000 


0.78336 


0.60000 


0.S6125 


0.70000 


0.93921 


O.SOOOO 


1.01720 


0.90000 


1.09520 


1.00000 


1.17320 


0. 10000 


0.5000S 


0.20000 


0.56605 


0.30000 


0.6329S 


0.40000 


0.70021 


0.50000 


0.76757 


0.60000 


0.83499 


0.70000 


0.90245 


0.80000 


0.96993 


0.90000 


1.03740 


1.00000 


1.10490 


0.10000 


0.531S9 


0.20000 


0.59036 


0.30000 


0.64946 


0.40000 


0.70875 


0.50000 


0.76813 


0.60000 


0.82756 



0.7599 


0.9005 


1.3845 


1.1390 


2.1813 


1.3762 


3.1372 


1.6120 


4.2318 


1.8455 


5.4501 


2.0760 


6.7009 


2.29S4 


8.0218 


2.5160 


9.4955 


2.7345 


0.2533 


0.6421 


0.5490 


0.8301 


0.9905 


1.02S7 


1.5579 


1.2275 


2.2462 


1.4262 


3.0514 


1.6245 


3.9528 


1.8211 


4.9258 


2.0147 


5.9872 


2.2067 


7.1377 


2.3975 


0.3S78 


0.7400 


0.4620 


0.8138 


0.7S05 


0.9736 


1.2014 


1.1408 


1.7177 


1.3100 


2.3256 


1.4S01 


3.0156 


1.6500 


3.7809 


1.8192 


4.6278 


1.9S82 


5.5291 


2.1552 


0.5459 


0.8400 


0.4757 


0.8496 


0.6791 


0.9597 


0.9953 


1.0967 


1.3928 


1.2406 


1.8649 


1.3873 
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o.soooo 

o.soooo 

0.80000 

o.soooo 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.6000 



0.70000 

0.80000 

0.90000 

1.00000 

0.10000 

0.20000 

0.30000 

0.40000 

0.50000 

0.60000 

0.70000 

0.80000 

0.90000 

1.00000 

0.10000 

0.20000 

0.30000 

0.30000 

0.40000 

0.50000 

0.60000 

0.70000 

O.SOOOO 

0.90000 

1.00000 



0.88701 

0.94648 

1.00600 

1.06540 

0.56816 

0.62062 

0.67352 

0.72655 

0.77964 

0. S 3276 

0.88590 

0.93905 

0.99221 

1.04540 

0.60743 

0.6549 S 

0.70284 

0.702 S 4 

0.75080 

0.79880 

0. S 4683 

0.89487 

0.94292 

0.99097 

1.14690 



2.4071 

3.0174 

3.6918 

4.4141 

0.7239 

0.6242 

0.6623 

0.8877 

1.1953 

1.5688 

2.0029 

2.4950 

3.0396 

3.6298 

0.9219 

0. S 218 

0.7427 

0.7427 

0.8598 

1.0SS6 

1.3838 

1.7344 

2.1352 

2.5 S 24 

2.0727 



1.5352 

1.6 S 38 

1.8324 

1.9801 

0.9400 

0.9400 

0.9825 

1.0862 

1.2059 

1.3317 

1.4606 

1.5913 

1.7226 

1.8540 

1.0400 

1.0400 

1.044 S 

1.044 S 

1.1062 

1.1999 

1.3055 

1.4169 

1.5314 

1.6476 

1.7147 
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E. SIMULATION PROGRAM FOR EXAMPLE 4-3. 



TITLE PERFORMANCE INDEX CALCULATION 
ARRAY A( 1) ,B(3) ,C(2) ,D(2) 

TABLE A(1)=0. 1084 , B( l-3)=90. 36, 1. 0,0. 0,D( 1)=1. 0 
PARAM Z=0. 0,P=0. 0,Y0=0. 0 
CONST K1=0. 454616, R=16. 296 
INITIAL 

C(1)=Z 

C(2)=P 

D(2)=P 

BETA=P 

ALPHA=Z 

DYNAMIC 

R=STEP( 0. ) 

DERIVATIVE 

SI=TRNFR( 0 , 2 , Y0 , A,B ,DEL) 

DEL=TRNFR( 1, 1,Y0,C,D,X1) 

SIE=SI -R 
X1=K1*SIE 
SIES=SIE*SIE 
DELS=DEL*DEL 
Q=R SIES+DELS 
J=INTGRL( YO ,Q) 

CONTROL FINTIM=15. ,DELT=0. 05 



PRINT 

END 


15. , 


ALPHA, BETA, J 


PARAM 

END 


Z=2. 


,P=0. 1 


PARAM 

END 


Z=2. 


,P=0. 3 


PARAM 

END 


Z=2. 


,P=0. 5 


PARAM 

END 


CM 

II 

CM 


,P=0. 7 


PARAM 

END 


Z=2. 


,P=0. 9 


PARAM 

END 


Z=2. 


,P=1. 0 


PARAM 

END 


II 

CM 


5 , P=0. 1 


PARAM 

END 


Z=4. 


5 , P=0. 3 


PARAM 

END 


II 

CM 


5 ,P=0. 5 


PARAM 

END 


Z=4. 


5 , P=0 . 7 


PARAM 

END 


CM 

II 


5 , P=0. 9 


PARAM 

END 


CM 

II 


5 ,P=1. 0 


PARAM 

END 


CM 

II 

^4 


5 , P=0. 1 


PARAM 

END 


CM 

II 


5 , P=0. 3 


PARAM 


Z=7. 


5 , P=0. 5 
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END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 



Z=7. 5,P=0. 7 
Z=7. 5 , P=0. 9 
Z=7. 5 , P=l. 0 
Z=9. ,P=0. 1 
Z=9. ,P=0. 3 
Z=9. ,P=0. 5 
Z=9. ,P=0. 7 
Z=9. , P=0. 9 
Z=9. ,P=1. 0 
Z=13. ,P=0. 1 
Z=13. ,P=0. 3 
Z=13. ,P=0. 5 
Z=13. ,P=0. 7 
Z=13. ,P=0. 9 
Z=13. ,P=1. 0 
Z=16. 5 ,P=0. 1 
Z=16. 5 ,P=0. 3 
Z=16. 5 ,P=0. 5 
Z=16. 5 ,P=0. 7 
Z=16. 5 , P=0. 9 
Z=16. 5 ,P=1. 0 
Z=19. ,P=0. 1 
Z=19. ,P=0. 3 
Z=19. ,P=0. 5 
Z=19. ,P=0. 7 
Z=19. ,P=0. 9 
Z=19. ,P=1. 0 
Z=23. 5 , P=0. 1 



END 






PARAM 

END 


Z=23. 


5 , P=0. 3 


PARAM 

END 


Z=23. 


5,P=0. 5 


PARAM 

END 


Z=23. 


5 , P=0. 7 


PARAM 

END 


CO 

CNI 

II 

NJ 


5,P=0. 9 


PARAM 

END 


CO 

eg 

II 

tS3 


5 ,P=1. 0 


PARAM 

END 


Z=26. 


,P=0. 1 


PARAM 

END 


Z=26. 


,P=0. 3 


PARAM 

END 


Z=26. 


,P=0. 5 


PARAM 

END 


Z=26. 


,P=0. 7 


PARAM 

END 


Z=26. 


, P=0. 9 


PARAM 

END 


Z=26. 


,P=1. 0 


PARAM 

END 


Z=29. 


5 , P=0. 1 


PARAM 

END 


Z=29. 


5 , P=0. 3 


PARAM 

END 


Z=29. 


5 , P=0. 5 


PARAM 

END 


Z=29. 


5 ,P=0. 7 


PARAM 

END 


Z=29. 


5 , P=0. 9 


PARAM 

END 


Z=29. 


5 ,P=1. 0 


PARAM 

END 


Z=35. 


,P=0. 1 


PARAM 

END 


Z=35. 


, P=0. 3 


PARAM 

END 


Z=35. 


,P=0. 5 


PARAM 

END 


Z=35. 


,P=0. 7 


PARAM 

END 


Z=35. 


,P=0. 9 


PARAM 

END 

STOP 


Z=35. 


,P=1. 0 
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F. DATA FILE FOR EXAMPLE 4-3. 

ALPHA BETA J 
0.0000 0.00000 15.000 

2.0000 0.10000 14.937 

2.0000 0.30000 15.050 

2.0000 0.50000 15.042 

2.0000 0.70000 15.020 

2.0000 0.90000 15.002 

2.0000 1.00000 14.994 

4.5000 0.10000 14.S99 

4.5000 0.30000 14.9S2 

4.5000 0.50000 14.982 

4.5000 0.70000 14.970 

4.5000 0.90000 14.959 

4.5000 1.00000 14.955 

7.5000 0.10000 14.S62 

7.5000 0.30000 14.912 

7.5000 0.50000 14.9 IS 

7.5000 0.90000 14.912 

7.5000 1.00000 14.911 

9.0000 0.10000 14.S46 

9.0000 0.30000 14. SS 1 

9.0000 0.50000 14.8S9 

9.0000 0.70000 14.S90 

9.0000 0.90000 14.890 

9.0000 1.00000 14.890 

13.0000 0.10000 14.812 

13.0000 0.30000 14.S09 

13.0000 0.50000 14.820 

13.0000 0.70000 14.829 

13.0000 0.90000 14.836 

13.0000 1.00000 14.S39 

16.5000 0.10000 14.791 
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16.5000 

16.5000 

16.5000 

16.5000 

16.5000 

19.0000 

19.0000 

19.000 

19.000 

19.000 

19.000 

23.500 

23.500 

23.500 

23.500 

23.500 

23.500 

26.000 

26.000 

26.000 

26.000 

26.000 

26.000 

29.500 

29.500 

29.500 

29.500 

29.500 

29.500 

35.000 

35.000 

35.000 

35.000 

35.000 

35.000 



0.30000 

0.50000 

0.70000 

0.90000 

1.00000 

0.10000 

0.30000 

0.50000 

0.70000 

0.90000 

1.00000 

0.10000 

0.30000 

0.50000 

0.70000 

0.90000 

1.00000 

0.10000 

0.30000 

0.50000 

0.70000 

0.90000 

1.00000 

0.10000 

0.30000 

0.50000 

0.70000 

0.90000 

1.00000 

1.00000 

0.30000 

0.50000 

0.70000 

0.90000 

1.00000 



14.759 
14.769 
14.782 
14.794 
14.799 
14.779 
14.730 
14.737 
14.752 
14.766 
14.773 
14.765 
14.686 
14.6 S 9 
14.704 
14.722 
14.730 

14.760 
14.667 
14.666 

14.6 5 1 
14.699 
14.708 
14.755 
14.646 
14.63 S 

14.652 
14.671 
14.680 
14.754 
14.621 
14.603 
14.614 
14.632 
14.642 



APPENDIX D. USER S MANUAL 



When we wish to plot the constant cost curves, the following described procedures 
should be performed as indicated. In this section, let us consider the specific problem 
in example 4-1. 

1. PLOTTING THE PARAMETER PLANE CURVES. 

The first step is to plot the parameter plane curves which are constant zeta curves 
and constant omega curves as shown in Figure 12 in Chapter 4. This curve permits us 
to choose the range of two parameters that are to be used in the analysis and design. 

2. GENERATING A DATA FILE. 



The second is to obtain a data file containing a , /? and performance index values. 
Simulation in D.S.L is used to evaluate the cost at a large number of points on the pa- 
rameter plane. The simulation results look as follows. 

1 DSL/VS RELEASE 1 I10D1F1CAT 1011 LEVEL 1 1985 



0 

o >>>> Dsi simulation mrur data <<<< 

TITLE PERTORHAIICE INDEX CALCULATIOII 

TADLE 2(I) = 1 00. ,P(1 -ATM. 0,11.0, 10. 0,0.0,DU) = 1.0 

FARAII A=0.O,D*0,0,Y0“0.0 

CUNTROl F 1 111 1 11= 1 5 » DEL 1 = 0.01 

mill 15.,ALr)IA,DEIA,lE5, 1 AET , 1AE 

EIID 

OXXX KLEIIGN = 2189, KPnillT = 227, AVAILABLE SPACE LEFT III COMMOIT/CURVAL/ * 

OXXX RK5 III! EI1RAT 1 Oil I1EIII0D USED xxx 

0 xxx DSI 0UIPU1 L1SI1IIG, GROUP 1 

+ , S TART I HO RUM 1 xxx 

PERFORITAIICE IIIDEX CALCULATION 

TIME ALPHA BEI A IES I AET IAE 

0.00000ET00 0 . 0O0 0 0E l 00 0 . 00000E* 00 O.OOOOOE+OO 0.00000E+00 0 . 00000E+00 

15.000 0.00000E»00 0.00000EI00 15.000 112.50 15.000 

0 

0 >>>> DSL SIMULATION 1IIPUT DATA <<<< 

PARAIT A = A. ,0 = 10. 

EIID 

0 xxx DSL OUIPUT L1ST1IIG, GROUP 1 

♦ , STARTING RUN 2 XXX 

PERFORMANCE IIIDEX CALCULATION 

TIME ALPHA DEI A IES 1 AET IAE 

0 . OOOOOE* 00 A. 0000 1 0.000 0.00000E+00 0.00000E+00 0.00000E + 00 

15.000 A. 0000 10.000 0.J2856 0.A7A62 0.6A1A0 



We have to edit the data file as shown in Appendix C to remove unnecessary parts. 
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3. LOADING A DATA FILE. 



The third step is to download this File From the mainframe computer to a floppy 
diskette. This data now on the floppy diskette can be copied on the IBM microcomputer 
under a designated Filename. In this example, the Filename is called pil.dat. Now, we 
are prepared to plot the families of performance index curves in the same range of pa- 
rameter values on the two parameter plane. 

4. PLOTTING THE COST CURVES. 



To run the performance index program, type in the executable file name, 
COSTFUN. The IBM microcomputers in the NPS Controls Lab are set up under the 
file management system 1DIR. To execute the program under the this system, enter the 
appropriate subdirectory’ by using the up/down arrow keys to position the cursor adja- 
cent to the subdirectory name, then depress the ENTER key. Use this same procedure 
to select COSTFUN.EXE and push ENTER to run the program. The preparation step 
for plotting the desired curves will appear. 



Enter proper values after each prompt as follows: 



1 . How many 


constant PI curves # 


= = = >6 




2 . Enter 


PI 


values 


, lwgt , ldig 


ff 


1===>0 . 27 


,1,2 


2 . Enter 


PI 


values 


, lwgt, ldig 




2===>0 . 2 , 


1,1 


2 . En te r 


PI 


values 


, lwgt , ldig 


# 


3 = = = >0 . 15 


,1,2 


2 . Enter 


PI 


values 


, lwgt , ldig 


« 


04 

o 

ii 

ii 

ii 


,1,2 


2 . En te r 


PI 


values 


, lwgt , ldig 


f» 


5===>0. 11 


,1,2 


2 . En te r 


PI 


values 


, lwgt , ldig 


# 


o 

/\ 

ti 

ii 

n 

CD 


1,1 


3 • En ter 


number of 


data file 


points ===>86 
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Next, adjust the plot size to the same parameter plane curves. 



*** This is the routine to Adjust the plot size *** 



Output plots are currently sized to fill the console screen. 

Would you like to adjust the plot size? 

Enter "y" or "n" ==> y 



A default plot factor size of 0.9 is used to fill the console screen. 
This plot size will be halved by entering a value of 0.45. 

It will be doubled with a factor entry of 1.8. 

Enter'factor size ==> 0.7 



The last step is to enter a data file name. 



*** This is the routine to LOAD data from a file *** 



File name should not exceed 8 characters 
File extension should not exceed 3 characters 



What is the file name (fn) and extension (ext)? 

Enter in form fn.ext (e.g. PARAPLAN . INP ) ===> pil.dat 
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After displaying the graph on the monitor, the program prompts you for a hardcopy 
of the graph. 

WOULD YOU LIKE HARDCOPY OF THESE PLOTS? (Y or N) 

y 



PRINTER/OUTPUT MENU 



PRINTER NO. 



1 j 

2 ! 

3 ! 

4 i 

5 1 

6 ! 

7 ! 

8 ! 

9 ! 

+++++++++++++++ 



PRINTER 



Epson FX-80 , All 
Epson FX-100, All 
Epson MX-100, All 
Epson RX-80, All 
Epson MX-80 k IBM Printer 
HP 7470A Graphics Plotter 
HP 7475A Graphics Plotter 
HP 758xB Series Plotters 
HP 2686A Laser Jet 
++++++++++++++++++++++++++++++ 



+ 



++++++++++++++ 



+ 



10 

11 

12 

99 



Graphics Monitor (default) 

HARDWARE Interface Menu 

Input PLOT88 Values for IOPORT and MODEL 
RETURN to LOAD Data 



Enter integer number for selection = = = > 5 

Reenter the data file name for the hard copy. Figure 13 is a performance index 
curve for constraint f E 2 dt. In like manner, the other two nonlinear criterion of per- 
formance index in terms of f\E\ t dt and f \ E\ di can also be plotted by running the 
same program. Figure 14 and Figure 15 are shown for desired constraint f\E\ t dt and 
f\E\ di , respectively. 
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